library(ggplot2)
library(dplyr)
##
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
library(tidyverse)
## ── Attaching packages ──────────────────────────────────────────────── tidyverse 1.3.0 ──
## ✓ tibble 3.0.0 ✓ purrr 0.3.4
## ✓ tidyr 1.0.2 ✓ stringr 1.4.0
## ✓ readr 1.3.1 ✓ forcats 0.5.0
## ── Conflicts ─────────────────────────────────────────────────── tidyverse_conflicts() ──
## x dplyr::filter() masks stats::filter()
## x dplyr::lag() masks stats::lag()
library(knitr)
library(gridExtra)
##
## Attaching package: 'gridExtra'
## The following object is masked from 'package:dplyr':
##
## combine
library(stringr)
library(purrr)
samples_sinus <- read_csv("samples_normal_sinus_rythmn.csv", skip = 1)
## Parsed with column specification:
## cols(
## `'seconds'` = col_double(),
## `'mV'` = col_double()
## )
colnames(samples_sinus) <- c("time", "signal")
S <- 5000
burn <- 5000
y <- samples_sinus$signal[55:210]
t <- samples_sinus$time[55:210]
beta_curr <- 1
alpha_curr <- rep(1, 12)
delta_curr<- seq(0, max(t), length.out = 17)
end_time <- max(t)
tao_curr <- 1
mean_curr <- rep(mean(y), length(y)) #mean_curr is a vector cause mean_star of y is a vector so I think r only lets me run it if it is a vector
dc_curr <- 0.5
miu_curr <- mean_curr - dc_curr
#storage
DC <- NULL
ALPHA <- rep(list(NULL), 12)
BETA <- NULL
DELTA <- rep(list(NULL), 17)
TAO <- NULL
MEAN <- rep(list(NULL), length(y))
for(s in 1:S) {
print(c("s: ", s))
#1 draw of Beta
beta_lst <- beta_one_samp(beta_curr, BETA, s)
beta_curr <- beta_lst[[1]]
BETA <- beta_lst[[2]]
#1 draw of Alphas (1-12)
for (k in 1:12) {
alpha_lst <- alpha_one_samp(k, alpha_curr[k], ALPHA[[k]], s)
alpha_curr[k] <- alpha_lst[[1]]
ALPHA[[k]] <- alpha_lst[[2]]
}
#1 draw of Deltas (1-17)
delta_lst <- delta_one_samp(delta_curr[1], 0, delta_curr[2], DELTA[[1]], s)
delta_curr[1] <- delta_lst[[1]]
DELTA[[1]] <- delta_lst[[2]]
for (k in 2:16) {
delta_lst <- delta_one_samp(delta_curr[k], delta_curr[k-1], delta_curr[k+1], DELTA[[k]], s)
delta_curr[k] <- delta_lst[[1]]
DELTA[[k]] <- delta_lst[[2]]
}
delta_lst <- delta_one_samp(delta_curr[17], delta_curr[[16]], end_time, DELTA[[17]], s)
delta_curr[17] <- delta_lst[[1]]
DELTA[[17]] <- delta_lst[[2]]
#1 draw of Tao
tao_lst <- tao_one_samp_cond(tao_curr, TAO, s)
tao_curr <- tao_lst[[1]]
TAO <- tao_lst[[2]]
#1 draw of DC
dc_lst <- dc_one_samp(dc_curr, DC, s)
dc_curr <- dc_lst[[1]]
DC <- dc_lst[[2]]
mean_curr <- dc_curr + gen_miu_star(t, alpha_curr = alpha_curr, beta_curr = beta_curr, delta_curr = delta_curr)
for (k in 1:length(y)) {
MEAN[[k]] <- c(MEAN[[k]], mean_curr[[k]])
}
#print(c("MEAN CURR: ", mean_curr))
}
## [1] "s: " "1"
## [1] "s: " "2"
## [1] "s: " "3"
## [1] "s: " "4"
## [1] "s: " "5"
## [1] "s: " "6"
## [1] "s: " "7"
## [1] "s: " "8"
## [1] "s: " "9"
## [1] "s: " "10"
## [1] "s: " "11"
## [1] "s: " "12"
## [1] "s: " "13"
## [1] "s: " "14"
## [1] "s: " "15"
## [1] "s: " "16"
## [1] "s: " "17"
## [1] "s: " "18"
## [1] "s: " "19"
## [1] "s: " "20"
## [1] "s: " "21"
## [1] "s: " "22"
## [1] "s: " "23"
## [1] "s: " "24"
## [1] "s: " "25"
## [1] "s: " "26"
## [1] "s: " "27"
## [1] "s: " "28"
## [1] "s: " "29"
## [1] "s: " "30"
## [1] "s: " "31"
## [1] "s: " "32"
## [1] "s: " "33"
## [1] "s: " "34"
## [1] "s: " "35"
## [1] "s: " "36"
## [1] "s: " "37"
## [1] "s: " "38"
## [1] "s: " "39"
## [1] "s: " "40"
## [1] "s: " "41"
## [1] "s: " "42"
## [1] "s: " "43"
## [1] "s: " "44"
## [1] "s: " "45"
## [1] "s: " "46"
## [1] "s: " "47"
## [1] "s: " "48"
## [1] "s: " "49"
## [1] "s: " "50"
## [1] "s: " "51"
## [1] "s: " "52"
## [1] "s: " "53"
## [1] "s: " "54"
## [1] "s: " "55"
## [1] "s: " "56"
## [1] "s: " "57"
## [1] "s: " "58"
## [1] "s: " "59"
## [1] "s: " "60"
## [1] "s: " "61"
## [1] "s: " "62"
## [1] "s: " "63"
## [1] "s: " "64"
## [1] "s: " "65"
## [1] "s: " "66"
## [1] "s: " "67"
## [1] "s: " "68"
## [1] "s: " "69"
## [1] "s: " "70"
## [1] "s: " "71"
## [1] "s: " "72"
## [1] "s: " "73"
## [1] "s: " "74"
## [1] "s: " "75"
## [1] "s: " "76"
## [1] "s: " "77"
## [1] "s: " "78"
## [1] "s: " "79"
## [1] "s: " "80"
## [1] "s: " "81"
## [1] "s: " "82"
## [1] "s: " "83"
## [1] "s: " "84"
## [1] "s: " "85"
## [1] "s: " "86"
## [1] "s: " "87"
## [1] "s: " "88"
## [1] "s: " "89"
## [1] "s: " "90"
## [1] "s: " "91"
## [1] "s: " "92"
## [1] "s: " "93"
## [1] "s: " "94"
## [1] "s: " "95"
## [1] "s: " "96"
## [1] "s: " "97"
## [1] "s: " "98"
## [1] "s: " "99"
## [1] "s: " "100"
## [1] "s: " "101"
## [1] "s: " "102"
## [1] "s: " "103"
## [1] "s: " "104"
## [1] "s: " "105"
## [1] "s: " "106"
## [1] "s: " "107"
## [1] "s: " "108"
## [1] "s: " "109"
## [1] "s: " "110"
## [1] "s: " "111"
## [1] "s: " "112"
## [1] "s: " "113"
## [1] "s: " "114"
## [1] "s: " "115"
## [1] "s: " "116"
## [1] "s: " "117"
## [1] "s: " "118"
## [1] "s: " "119"
## [1] "s: " "120"
## [1] "s: " "121"
## [1] "s: " "122"
## [1] "s: " "123"
## [1] "s: " "124"
## [1] "s: " "125"
## [1] "s: " "126"
## [1] "s: " "127"
## [1] "s: " "128"
## [1] "s: " "129"
## [1] "s: " "130"
## [1] "s: " "131"
## [1] "s: " "132"
## [1] "s: " "133"
## [1] "s: " "134"
## [1] "s: " "135"
## [1] "s: " "136"
## [1] "s: " "137"
## [1] "s: " "138"
## [1] "s: " "139"
## [1] "s: " "140"
## [1] "s: " "141"
## [1] "s: " "142"
## [1] "s: " "143"
## [1] "s: " "144"
## [1] "s: " "145"
## [1] "s: " "146"
## [1] "s: " "147"
## [1] "s: " "148"
## [1] "s: " "149"
## [1] "s: " "150"
## [1] "s: " "151"
## [1] "s: " "152"
## [1] "s: " "153"
## [1] "s: " "154"
## [1] "s: " "155"
## [1] "s: " "156"
## [1] "s: " "157"
## [1] "s: " "158"
## [1] "s: " "159"
## [1] "s: " "160"
## [1] "s: " "161"
## [1] "s: " "162"
## [1] "s: " "163"
## [1] "s: " "164"
## [1] "s: " "165"
## [1] "s: " "166"
## [1] "s: " "167"
## [1] "s: " "168"
## [1] "s: " "169"
## [1] "s: " "170"
## [1] "s: " "171"
## [1] "s: " "172"
## [1] "s: " "173"
## [1] "s: " "174"
## [1] "s: " "175"
## [1] "s: " "176"
## [1] "s: " "177"
## [1] "s: " "178"
## [1] "s: " "179"
## [1] "s: " "180"
## [1] "s: " "181"
## [1] "s: " "182"
## [1] "s: " "183"
## [1] "s: " "184"
## [1] "s: " "185"
## [1] "s: " "186"
## [1] "s: " "187"
## [1] "s: " "188"
## [1] "s: " "189"
## [1] "s: " "190"
## [1] "s: " "191"
## [1] "s: " "192"
## [1] "s: " "193"
## [1] "s: " "194"
## [1] "s: " "195"
## [1] "s: " "196"
## [1] "s: " "197"
## [1] "s: " "198"
## [1] "s: " "199"
## [1] "s: " "200"
## [1] "s: " "201"
## [1] "s: " "202"
## [1] "s: " "203"
## [1] "s: " "204"
## [1] "s: " "205"
## [1] "s: " "206"
## [1] "s: " "207"
## [1] "s: " "208"
## [1] "s: " "209"
## [1] "s: " "210"
## [1] "s: " "211"
## [1] "s: " "212"
## [1] "s: " "213"
## [1] "s: " "214"
## [1] "s: " "215"
## [1] "s: " "216"
## [1] "s: " "217"
## [1] "s: " "218"
## [1] "s: " "219"
## [1] "s: " "220"
## [1] "s: " "221"
## [1] "s: " "222"
## [1] "s: " "223"
## [1] "s: " "224"
## [1] "s: " "225"
## [1] "s: " "226"
## [1] "s: " "227"
## [1] "s: " "228"
## [1] "s: " "229"
## [1] "s: " "230"
## [1] "s: " "231"
## [1] "s: " "232"
## [1] "s: " "233"
## [1] "s: " "234"
## [1] "s: " "235"
## [1] "s: " "236"
## [1] "s: " "237"
## [1] "s: " "238"
## [1] "s: " "239"
## [1] "s: " "240"
## [1] "s: " "241"
## [1] "s: " "242"
## [1] "s: " "243"
## [1] "s: " "244"
## [1] "s: " "245"
## [1] "s: " "246"
## [1] "s: " "247"
## [1] "s: " "248"
## [1] "s: " "249"
## [1] "s: " "250"
## [1] "s: " "251"
## [1] "s: " "252"
## [1] "s: " "253"
## [1] "s: " "254"
## [1] "s: " "255"
## [1] "s: " "256"
## [1] "s: " "257"
## [1] "s: " "258"
## [1] "s: " "259"
## [1] "s: " "260"
## [1] "s: " "261"
## [1] "s: " "262"
## [1] "s: " "263"
## [1] "s: " "264"
## [1] "s: " "265"
## [1] "s: " "266"
## [1] "s: " "267"
## [1] "s: " "268"
## [1] "s: " "269"
## [1] "s: " "270"
## [1] "s: " "271"
## [1] "s: " "272"
## [1] "s: " "273"
## [1] "s: " "274"
## [1] "s: " "275"
## [1] "s: " "276"
## [1] "s: " "277"
## [1] "s: " "278"
## [1] "s: " "279"
## [1] "s: " "280"
## [1] "s: " "281"
## [1] "s: " "282"
## [1] "s: " "283"
## [1] "s: " "284"
## [1] "s: " "285"
## [1] "s: " "286"
## [1] "s: " "287"
## [1] "s: " "288"
## [1] "s: " "289"
## [1] "s: " "290"
## [1] "s: " "291"
## [1] "s: " "292"
## [1] "s: " "293"
## [1] "s: " "294"
## [1] "s: " "295"
## [1] "s: " "296"
## [1] "s: " "297"
## [1] "s: " "298"
## [1] "s: " "299"
## [1] "s: " "300"
## [1] "s: " "301"
## [1] "s: " "302"
## [1] "s: " "303"
## [1] "s: " "304"
## [1] "s: " "305"
## [1] "s: " "306"
## [1] "s: " "307"
## [1] "s: " "308"
## [1] "s: " "309"
## [1] "s: " "310"
## [1] "s: " "311"
## [1] "s: " "312"
## [1] "s: " "313"
## [1] "s: " "314"
## [1] "s: " "315"
## [1] "s: " "316"
## [1] "s: " "317"
## [1] "s: " "318"
## [1] "s: " "319"
## [1] "s: " "320"
## [1] "s: " "321"
## [1] "s: " "322"
## [1] "s: " "323"
## [1] "s: " "324"
## [1] "s: " "325"
## [1] "s: " "326"
## [1] "s: " "327"
## [1] "s: " "328"
## [1] "s: " "329"
## [1] "s: " "330"
## [1] "s: " "331"
## [1] "s: " "332"
## [1] "s: " "333"
## [1] "s: " "334"
## [1] "s: " "335"
## [1] "s: " "336"
## [1] "s: " "337"
## [1] "s: " "338"
## [1] "s: " "339"
## [1] "s: " "340"
## [1] "s: " "341"
## [1] "s: " "342"
## [1] "s: " "343"
## [1] "s: " "344"
## [1] "s: " "345"
## [1] "s: " "346"
## [1] "s: " "347"
## [1] "s: " "348"
## [1] "s: " "349"
## [1] "s: " "350"
## [1] "s: " "351"
## [1] "s: " "352"
## [1] "s: " "353"
## [1] "s: " "354"
## [1] "s: " "355"
## [1] "s: " "356"
## [1] "s: " "357"
## [1] "s: " "358"
## [1] "s: " "359"
## [1] "s: " "360"
## [1] "s: " "361"
## [1] "s: " "362"
## [1] "s: " "363"
## [1] "s: " "364"
## [1] "s: " "365"
## [1] "s: " "366"
## [1] "s: " "367"
## [1] "s: " "368"
## [1] "s: " "369"
## [1] "s: " "370"
## [1] "s: " "371"
## [1] "s: " "372"
## [1] "s: " "373"
## [1] "s: " "374"
## [1] "s: " "375"
## [1] "s: " "376"
## [1] "s: " "377"
## [1] "s: " "378"
## [1] "s: " "379"
## [1] "s: " "380"
## [1] "s: " "381"
## [1] "s: " "382"
## [1] "s: " "383"
## [1] "s: " "384"
## [1] "s: " "385"
## [1] "s: " "386"
## [1] "s: " "387"
## [1] "s: " "388"
## [1] "s: " "389"
## [1] "s: " "390"
## [1] "s: " "391"
## [1] "s: " "392"
## [1] "s: " "393"
## [1] "s: " "394"
## [1] "s: " "395"
## [1] "s: " "396"
## [1] "s: " "397"
## [1] "s: " "398"
## [1] "s: " "399"
## [1] "s: " "400"
## [1] "s: " "401"
## [1] "s: " "402"
## [1] "s: " "403"
## [1] "s: " "404"
## [1] "s: " "405"
## [1] "s: " "406"
## [1] "s: " "407"
## [1] "s: " "408"
## [1] "s: " "409"
## [1] "s: " "410"
## [1] "s: " "411"
## [1] "s: " "412"
## [1] "s: " "413"
## [1] "s: " "414"
## [1] "s: " "415"
## [1] "s: " "416"
## [1] "s: " "417"
## [1] "s: " "418"
## [1] "s: " "419"
## [1] "s: " "420"
## [1] "s: " "421"
## [1] "s: " "422"
## [1] "s: " "423"
## [1] "s: " "424"
## [1] "s: " "425"
## [1] "s: " "426"
## [1] "s: " "427"
## [1] "s: " "428"
## [1] "s: " "429"
## [1] "s: " "430"
## [1] "s: " "431"
## [1] "s: " "432"
## [1] "s: " "433"
## [1] "s: " "434"
## [1] "s: " "435"
## [1] "s: " "436"
## [1] "s: " "437"
## [1] "s: " "438"
## [1] "s: " "439"
## [1] "s: " "440"
## [1] "s: " "441"
## [1] "s: " "442"
## [1] "s: " "443"
## [1] "s: " "444"
## [1] "s: " "445"
## [1] "s: " "446"
## [1] "s: " "447"
## [1] "s: " "448"
## [1] "s: " "449"
## [1] "s: " "450"
## [1] "s: " "451"
## [1] "s: " "452"
## [1] "s: " "453"
## [1] "s: " "454"
## [1] "s: " "455"
## [1] "s: " "456"
## [1] "s: " "457"
## [1] "s: " "458"
## [1] "s: " "459"
## [1] "s: " "460"
## [1] "s: " "461"
## [1] "s: " "462"
## [1] "s: " "463"
## [1] "s: " "464"
## [1] "s: " "465"
## [1] "s: " "466"
## [1] "s: " "467"
## [1] "s: " "468"
## [1] "s: " "469"
## [1] "s: " "470"
## [1] "s: " "471"
## [1] "s: " "472"
## [1] "s: " "473"
## [1] "s: " "474"
## [1] "s: " "475"
## [1] "s: " "476"
## [1] "s: " "477"
## [1] "s: " "478"
## [1] "s: " "479"
## [1] "s: " "480"
## [1] "s: " "481"
## [1] "s: " "482"
## [1] "s: " "483"
## [1] "s: " "484"
## [1] "s: " "485"
## [1] "s: " "486"
## [1] "s: " "487"
## [1] "s: " "488"
## [1] "s: " "489"
## [1] "s: " "490"
## [1] "s: " "491"
## [1] "s: " "492"
## [1] "s: " "493"
## [1] "s: " "494"
## [1] "s: " "495"
## [1] "s: " "496"
## [1] "s: " "497"
## [1] "s: " "498"
## [1] "s: " "499"
## [1] "s: " "500"
## [1] "s: " "501"
## [1] "s: " "502"
## [1] "s: " "503"
## [1] "s: " "504"
## [1] "s: " "505"
## [1] "s: " "506"
## [1] "s: " "507"
## [1] "s: " "508"
## [1] "s: " "509"
## [1] "s: " "510"
## [1] "s: " "511"
## [1] "s: " "512"
## [1] "s: " "513"
## [1] "s: " "514"
## [1] "s: " "515"
## [1] "s: " "516"
## [1] "s: " "517"
## [1] "s: " "518"
## [1] "s: " "519"
## [1] "s: " "520"
## [1] "s: " "521"
## [1] "s: " "522"
## [1] "s: " "523"
## [1] "s: " "524"
## [1] "s: " "525"
## [1] "s: " "526"
## [1] "s: " "527"
## [1] "s: " "528"
## [1] "s: " "529"
## [1] "s: " "530"
## [1] "s: " "531"
## [1] "s: " "532"
## [1] "s: " "533"
## [1] "s: " "534"
## [1] "s: " "535"
## [1] "s: " "536"
## [1] "s: " "537"
## [1] "s: " "538"
## [1] "s: " "539"
## [1] "s: " "540"
## [1] "s: " "541"
## [1] "s: " "542"
## [1] "s: " "543"
## [1] "s: " "544"
## [1] "s: " "545"
## [1] "s: " "546"
## [1] "s: " "547"
## [1] "s: " "548"
## [1] "s: " "549"
## [1] "s: " "550"
## [1] "s: " "551"
## [1] "s: " "552"
## [1] "s: " "553"
## [1] "s: " "554"
## [1] "s: " "555"
## [1] "s: " "556"
## [1] "s: " "557"
## [1] "s: " "558"
## [1] "s: " "559"
## [1] "s: " "560"
## [1] "s: " "561"
## [1] "s: " "562"
## [1] "s: " "563"
## [1] "s: " "564"
## [1] "s: " "565"
## [1] "s: " "566"
## [1] "s: " "567"
## [1] "s: " "568"
## [1] "s: " "569"
## [1] "s: " "570"
## [1] "s: " "571"
## [1] "s: " "572"
## [1] "s: " "573"
## [1] "s: " "574"
## [1] "s: " "575"
## [1] "s: " "576"
## [1] "s: " "577"
## [1] "s: " "578"
## [1] "s: " "579"
## [1] "s: " "580"
## [1] "s: " "581"
## [1] "s: " "582"
## [1] "s: " "583"
## [1] "s: " "584"
## [1] "s: " "585"
## [1] "s: " "586"
## [1] "s: " "587"
## [1] "s: " "588"
## [1] "s: " "589"
## [1] "s: " "590"
## [1] "s: " "591"
## [1] "s: " "592"
## [1] "s: " "593"
## [1] "s: " "594"
## [1] "s: " "595"
## [1] "s: " "596"
## [1] "s: " "597"
## [1] "s: " "598"
## [1] "s: " "599"
## [1] "s: " "600"
## [1] "s: " "601"
## [1] "s: " "602"
## [1] "s: " "603"
## [1] "s: " "604"
## [1] "s: " "605"
## [1] "s: " "606"
## [1] "s: " "607"
## [1] "s: " "608"
## [1] "s: " "609"
## [1] "s: " "610"
## [1] "s: " "611"
## [1] "s: " "612"
## [1] "s: " "613"
## [1] "s: " "614"
## [1] "s: " "615"
## [1] "s: " "616"
## [1] "s: " "617"
## [1] "s: " "618"
## [1] "s: " "619"
## [1] "s: " "620"
## [1] "s: " "621"
## [1] "s: " "622"
## [1] "s: " "623"
## [1] "s: " "624"
## [1] "s: " "625"
## [1] "s: " "626"
## [1] "s: " "627"
## [1] "s: " "628"
## [1] "s: " "629"
## [1] "s: " "630"
## [1] "s: " "631"
## [1] "s: " "632"
## [1] "s: " "633"
## [1] "s: " "634"
## [1] "s: " "635"
## [1] "s: " "636"
## [1] "s: " "637"
## [1] "s: " "638"
## [1] "s: " "639"
## [1] "s: " "640"
## [1] "s: " "641"
## [1] "s: " "642"
## [1] "s: " "643"
## [1] "s: " "644"
## [1] "s: " "645"
## [1] "s: " "646"
## [1] "s: " "647"
## [1] "s: " "648"
## [1] "s: " "649"
## [1] "s: " "650"
## [1] "s: " "651"
## [1] "s: " "652"
## [1] "s: " "653"
## [1] "s: " "654"
## [1] "s: " "655"
## [1] "s: " "656"
## [1] "s: " "657"
## [1] "s: " "658"
## [1] "s: " "659"
## [1] "s: " "660"
## [1] "s: " "661"
## [1] "s: " "662"
## [1] "s: " "663"
## [1] "s: " "664"
## [1] "s: " "665"
## [1] "s: " "666"
## [1] "s: " "667"
## [1] "s: " "668"
## [1] "s: " "669"
## [1] "s: " "670"
## [1] "s: " "671"
## [1] "s: " "672"
## [1] "s: " "673"
## [1] "s: " "674"
## [1] "s: " "675"
## [1] "s: " "676"
## [1] "s: " "677"
## [1] "s: " "678"
## [1] "s: " "679"
## [1] "s: " "680"
## [1] "s: " "681"
## [1] "s: " "682"
## [1] "s: " "683"
## [1] "s: " "684"
## [1] "s: " "685"
## [1] "s: " "686"
## [1] "s: " "687"
## [1] "s: " "688"
## [1] "s: " "689"
## [1] "s: " "690"
## [1] "s: " "691"
## [1] "s: " "692"
## [1] "s: " "693"
## [1] "s: " "694"
## [1] "s: " "695"
## [1] "s: " "696"
## [1] "s: " "697"
## [1] "s: " "698"
## [1] "s: " "699"
## [1] "s: " "700"
## [1] "s: " "701"
## [1] "s: " "702"
## [1] "s: " "703"
## [1] "s: " "704"
## [1] "s: " "705"
## [1] "s: " "706"
## [1] "s: " "707"
## [1] "s: " "708"
## [1] "s: " "709"
## [1] "s: " "710"
## [1] "s: " "711"
## [1] "s: " "712"
## [1] "s: " "713"
## [1] "s: " "714"
## [1] "s: " "715"
## [1] "s: " "716"
## [1] "s: " "717"
## [1] "s: " "718"
## [1] "s: " "719"
## [1] "s: " "720"
## [1] "s: " "721"
## [1] "s: " "722"
## [1] "s: " "723"
## [1] "s: " "724"
## [1] "s: " "725"
## [1] "s: " "726"
## [1] "s: " "727"
## [1] "s: " "728"
## [1] "s: " "729"
## [1] "s: " "730"
## [1] "s: " "731"
## [1] "s: " "732"
## [1] "s: " "733"
## [1] "s: " "734"
## [1] "s: " "735"
## [1] "s: " "736"
## [1] "s: " "737"
## [1] "s: " "738"
## [1] "s: " "739"
## [1] "s: " "740"
## [1] "s: " "741"
## [1] "s: " "742"
## [1] "s: " "743"
## [1] "s: " "744"
## [1] "s: " "745"
## [1] "s: " "746"
## [1] "s: " "747"
## [1] "s: " "748"
## [1] "s: " "749"
## [1] "s: " "750"
## [1] "s: " "751"
## [1] "s: " "752"
## [1] "s: " "753"
## [1] "s: " "754"
## [1] "s: " "755"
## [1] "s: " "756"
## [1] "s: " "757"
## [1] "s: " "758"
## [1] "s: " "759"
## [1] "s: " "760"
## [1] "s: " "761"
## [1] "s: " "762"
## [1] "s: " "763"
## [1] "s: " "764"
## [1] "s: " "765"
## [1] "s: " "766"
## [1] "s: " "767"
## [1] "s: " "768"
## [1] "s: " "769"
## [1] "s: " "770"
## [1] "s: " "771"
## [1] "s: " "772"
## [1] "s: " "773"
## [1] "s: " "774"
## [1] "s: " "775"
## [1] "s: " "776"
## [1] "s: " "777"
## [1] "s: " "778"
## [1] "s: " "779"
## [1] "s: " "780"
## [1] "s: " "781"
## [1] "s: " "782"
## [1] "s: " "783"
## [1] "s: " "784"
## [1] "s: " "785"
## [1] "s: " "786"
## [1] "s: " "787"
## [1] "s: " "788"
## [1] "s: " "789"
## [1] "s: " "790"
## [1] "s: " "791"
## [1] "s: " "792"
## [1] "s: " "793"
## [1] "s: " "794"
## [1] "s: " "795"
## [1] "s: " "796"
## [1] "s: " "797"
## [1] "s: " "798"
## [1] "s: " "799"
## [1] "s: " "800"
## [1] "s: " "801"
## [1] "s: " "802"
## [1] "s: " "803"
## [1] "s: " "804"
## [1] "s: " "805"
## [1] "s: " "806"
## [1] "s: " "807"
## [1] "s: " "808"
## [1] "s: " "809"
## [1] "s: " "810"
## [1] "s: " "811"
## [1] "s: " "812"
## [1] "s: " "813"
## [1] "s: " "814"
## [1] "s: " "815"
## [1] "s: " "816"
## [1] "s: " "817"
## [1] "s: " "818"
## [1] "s: " "819"
## [1] "s: " "820"
## [1] "s: " "821"
## [1] "s: " "822"
## [1] "s: " "823"
## [1] "s: " "824"
## [1] "s: " "825"
## [1] "s: " "826"
## [1] "s: " "827"
## [1] "s: " "828"
## [1] "s: " "829"
## [1] "s: " "830"
## [1] "s: " "831"
## [1] "s: " "832"
## [1] "s: " "833"
## [1] "s: " "834"
## [1] "s: " "835"
## [1] "s: " "836"
## [1] "s: " "837"
## [1] "s: " "838"
## [1] "s: " "839"
## [1] "s: " "840"
## [1] "s: " "841"
## [1] "s: " "842"
## [1] "s: " "843"
## [1] "s: " "844"
## [1] "s: " "845"
## [1] "s: " "846"
## [1] "s: " "847"
## [1] "s: " "848"
## [1] "s: " "849"
## [1] "s: " "850"
## [1] "s: " "851"
## [1] "s: " "852"
## [1] "s: " "853"
## [1] "s: " "854"
## [1] "s: " "855"
## [1] "s: " "856"
## [1] "s: " "857"
## [1] "s: " "858"
## [1] "s: " "859"
## [1] "s: " "860"
## [1] "s: " "861"
## [1] "s: " "862"
## [1] "s: " "863"
## [1] "s: " "864"
## [1] "s: " "865"
## [1] "s: " "866"
## [1] "s: " "867"
## [1] "s: " "868"
## [1] "s: " "869"
## [1] "s: " "870"
## [1] "s: " "871"
## [1] "s: " "872"
## [1] "s: " "873"
## [1] "s: " "874"
## [1] "s: " "875"
## [1] "s: " "876"
## [1] "s: " "877"
## [1] "s: " "878"
## [1] "s: " "879"
## [1] "s: " "880"
## [1] "s: " "881"
## [1] "s: " "882"
## [1] "s: " "883"
## [1] "s: " "884"
## [1] "s: " "885"
## [1] "s: " "886"
## [1] "s: " "887"
## [1] "s: " "888"
## [1] "s: " "889"
## [1] "s: " "890"
## [1] "s: " "891"
## [1] "s: " "892"
## [1] "s: " "893"
## [1] "s: " "894"
## [1] "s: " "895"
## [1] "s: " "896"
## [1] "s: " "897"
## [1] "s: " "898"
## [1] "s: " "899"
## [1] "s: " "900"
## [1] "s: " "901"
## [1] "s: " "902"
## [1] "s: " "903"
## [1] "s: " "904"
## [1] "s: " "905"
## [1] "s: " "906"
## [1] "s: " "907"
## [1] "s: " "908"
## [1] "s: " "909"
## [1] "s: " "910"
## [1] "s: " "911"
## [1] "s: " "912"
## [1] "s: " "913"
## [1] "s: " "914"
## [1] "s: " "915"
## [1] "s: " "916"
## [1] "s: " "917"
## [1] "s: " "918"
## [1] "s: " "919"
## [1] "s: " "920"
## [1] "s: " "921"
## [1] "s: " "922"
## [1] "s: " "923"
## [1] "s: " "924"
## [1] "s: " "925"
## [1] "s: " "926"
## [1] "s: " "927"
## [1] "s: " "928"
## [1] "s: " "929"
## [1] "s: " "930"
## [1] "s: " "931"
## [1] "s: " "932"
## [1] "s: " "933"
## [1] "s: " "934"
## [1] "s: " "935"
## [1] "s: " "936"
## [1] "s: " "937"
## [1] "s: " "938"
## [1] "s: " "939"
## [1] "s: " "940"
## [1] "s: " "941"
## [1] "s: " "942"
## [1] "s: " "943"
## [1] "s: " "944"
## [1] "s: " "945"
## [1] "s: " "946"
## [1] "s: " "947"
## [1] "s: " "948"
## [1] "s: " "949"
## [1] "s: " "950"
## [1] "s: " "951"
## [1] "s: " "952"
## [1] "s: " "953"
## [1] "s: " "954"
## [1] "s: " "955"
## [1] "s: " "956"
## [1] "s: " "957"
## [1] "s: " "958"
## [1] "s: " "959"
## [1] "s: " "960"
## [1] "s: " "961"
## [1] "s: " "962"
## [1] "s: " "963"
## [1] "s: " "964"
## [1] "s: " "965"
## [1] "s: " "966"
## [1] "s: " "967"
## [1] "s: " "968"
## [1] "s: " "969"
## [1] "s: " "970"
## [1] "s: " "971"
## [1] "s: " "972"
## [1] "s: " "973"
## [1] "s: " "974"
## [1] "s: " "975"
## [1] "s: " "976"
## [1] "s: " "977"
## [1] "s: " "978"
## [1] "s: " "979"
## [1] "s: " "980"
## [1] "s: " "981"
## [1] "s: " "982"
## [1] "s: " "983"
## [1] "s: " "984"
## [1] "s: " "985"
## [1] "s: " "986"
## [1] "s: " "987"
## [1] "s: " "988"
## [1] "s: " "989"
## [1] "s: " "990"
## [1] "s: " "991"
## [1] "s: " "992"
## [1] "s: " "993"
## [1] "s: " "994"
## [1] "s: " "995"
## [1] "s: " "996"
## [1] "s: " "997"
## [1] "s: " "998"
## [1] "s: " "999"
## [1] "s: " "1000"
## [1] "s: " "1001"
## [1] "s: " "1002"
## [1] "s: " "1003"
## [1] "s: " "1004"
## [1] "s: " "1005"
## [1] "s: " "1006"
## [1] "s: " "1007"
## [1] "s: " "1008"
## [1] "s: " "1009"
## [1] "s: " "1010"
## [1] "s: " "1011"
## [1] "s: " "1012"
## [1] "s: " "1013"
## [1] "s: " "1014"
## [1] "s: " "1015"
## [1] "s: " "1016"
## [1] "s: " "1017"
## [1] "s: " "1018"
## [1] "s: " "1019"
## [1] "s: " "1020"
## [1] "s: " "1021"
## [1] "s: " "1022"
## [1] "s: " "1023"
## [1] "s: " "1024"
## [1] "s: " "1025"
## [1] "s: " "1026"
## [1] "s: " "1027"
## [1] "s: " "1028"
## [1] "s: " "1029"
## [1] "s: " "1030"
## [1] "s: " "1031"
## [1] "s: " "1032"
## [1] "s: " "1033"
## [1] "s: " "1034"
## [1] "s: " "1035"
## [1] "s: " "1036"
## [1] "s: " "1037"
## [1] "s: " "1038"
## [1] "s: " "1039"
## [1] "s: " "1040"
## [1] "s: " "1041"
## [1] "s: " "1042"
## [1] "s: " "1043"
## [1] "s: " "1044"
## [1] "s: " "1045"
## [1] "s: " "1046"
## [1] "s: " "1047"
## [1] "s: " "1048"
## [1] "s: " "1049"
## [1] "s: " "1050"
## [1] "s: " "1051"
## [1] "s: " "1052"
## [1] "s: " "1053"
## [1] "s: " "1054"
## [1] "s: " "1055"
## [1] "s: " "1056"
## [1] "s: " "1057"
## [1] "s: " "1058"
## [1] "s: " "1059"
## [1] "s: " "1060"
## [1] "s: " "1061"
## [1] "s: " "1062"
## [1] "s: " "1063"
## [1] "s: " "1064"
## [1] "s: " "1065"
## [1] "s: " "1066"
## [1] "s: " "1067"
## [1] "s: " "1068"
## [1] "s: " "1069"
## [1] "s: " "1070"
## [1] "s: " "1071"
## [1] "s: " "1072"
## [1] "s: " "1073"
## [1] "s: " "1074"
## [1] "s: " "1075"
## [1] "s: " "1076"
## [1] "s: " "1077"
## [1] "s: " "1078"
## [1] "s: " "1079"
## [1] "s: " "1080"
## [1] "s: " "1081"
## [1] "s: " "1082"
## [1] "s: " "1083"
## [1] "s: " "1084"
## [1] "s: " "1085"
## [1] "s: " "1086"
## [1] "s: " "1087"
## [1] "s: " "1088"
## [1] "s: " "1089"
## [1] "s: " "1090"
## [1] "s: " "1091"
## [1] "s: " "1092"
## [1] "s: " "1093"
## [1] "s: " "1094"
## [1] "s: " "1095"
## [1] "s: " "1096"
## [1] "s: " "1097"
## [1] "s: " "1098"
## [1] "s: " "1099"
## [1] "s: " "1100"
## [1] "s: " "1101"
## [1] "s: " "1102"
## [1] "s: " "1103"
## [1] "s: " "1104"
## [1] "s: " "1105"
## [1] "s: " "1106"
## [1] "s: " "1107"
## [1] "s: " "1108"
## [1] "s: " "1109"
## [1] "s: " "1110"
## [1] "s: " "1111"
## [1] "s: " "1112"
## [1] "s: " "1113"
## [1] "s: " "1114"
## [1] "s: " "1115"
## [1] "s: " "1116"
## [1] "s: " "1117"
## [1] "s: " "1118"
## [1] "s: " "1119"
## [1] "s: " "1120"
## [1] "s: " "1121"
## [1] "s: " "1122"
## [1] "s: " "1123"
## [1] "s: " "1124"
## [1] "s: " "1125"
## [1] "s: " "1126"
## [1] "s: " "1127"
## [1] "s: " "1128"
## [1] "s: " "1129"
## [1] "s: " "1130"
## [1] "s: " "1131"
## [1] "s: " "1132"
## [1] "s: " "1133"
## [1] "s: " "1134"
## [1] "s: " "1135"
## [1] "s: " "1136"
## [1] "s: " "1137"
## [1] "s: " "1138"
## [1] "s: " "1139"
## [1] "s: " "1140"
## [1] "s: " "1141"
## [1] "s: " "1142"
## [1] "s: " "1143"
## [1] "s: " "1144"
## [1] "s: " "1145"
## [1] "s: " "1146"
## [1] "s: " "1147"
## [1] "s: " "1148"
## [1] "s: " "1149"
## [1] "s: " "1150"
## [1] "s: " "1151"
## [1] "s: " "1152"
## [1] "s: " "1153"
## [1] "s: " "1154"
## [1] "s: " "1155"
## [1] "s: " "1156"
## [1] "s: " "1157"
## [1] "s: " "1158"
## [1] "s: " "1159"
## [1] "s: " "1160"
## [1] "s: " "1161"
## [1] "s: " "1162"
## [1] "s: " "1163"
## [1] "s: " "1164"
## [1] "s: " "1165"
## [1] "s: " "1166"
## [1] "s: " "1167"
## [1] "s: " "1168"
## [1] "s: " "1169"
## [1] "s: " "1170"
## [1] "s: " "1171"
## [1] "s: " "1172"
## [1] "s: " "1173"
## [1] "s: " "1174"
## [1] "s: " "1175"
## [1] "s: " "1176"
## [1] "s: " "1177"
## [1] "s: " "1178"
## [1] "s: " "1179"
## [1] "s: " "1180"
## [1] "s: " "1181"
## [1] "s: " "1182"
## [1] "s: " "1183"
## [1] "s: " "1184"
## [1] "s: " "1185"
## [1] "s: " "1186"
## [1] "s: " "1187"
## [1] "s: " "1188"
## [1] "s: " "1189"
## [1] "s: " "1190"
## [1] "s: " "1191"
## [1] "s: " "1192"
## [1] "s: " "1193"
## [1] "s: " "1194"
## [1] "s: " "1195"
## [1] "s: " "1196"
## [1] "s: " "1197"
## [1] "s: " "1198"
## [1] "s: " "1199"
## [1] "s: " "1200"
## [1] "s: " "1201"
## [1] "s: " "1202"
## [1] "s: " "1203"
## [1] "s: " "1204"
## [1] "s: " "1205"
## [1] "s: " "1206"
## [1] "s: " "1207"
## [1] "s: " "1208"
## [1] "s: " "1209"
## [1] "s: " "1210"
## [1] "s: " "1211"
## [1] "s: " "1212"
## [1] "s: " "1213"
## [1] "s: " "1214"
## [1] "s: " "1215"
## [1] "s: " "1216"
## [1] "s: " "1217"
## [1] "s: " "1218"
## [1] "s: " "1219"
## [1] "s: " "1220"
## [1] "s: " "1221"
## [1] "s: " "1222"
## [1] "s: " "1223"
## [1] "s: " "1224"
## [1] "s: " "1225"
## [1] "s: " "1226"
## [1] "s: " "1227"
## [1] "s: " "1228"
## [1] "s: " "1229"
## [1] "s: " "1230"
## [1] "s: " "1231"
## [1] "s: " "1232"
## [1] "s: " "1233"
## [1] "s: " "1234"
## [1] "s: " "1235"
## [1] "s: " "1236"
## [1] "s: " "1237"
## [1] "s: " "1238"
## [1] "s: " "1239"
## [1] "s: " "1240"
## [1] "s: " "1241"
## [1] "s: " "1242"
## [1] "s: " "1243"
## [1] "s: " "1244"
## [1] "s: " "1245"
## [1] "s: " "1246"
## [1] "s: " "1247"
## [1] "s: " "1248"
## [1] "s: " "1249"
## [1] "s: " "1250"
## [1] "s: " "1251"
## [1] "s: " "1252"
## [1] "s: " "1253"
## [1] "s: " "1254"
## [1] "s: " "1255"
## [1] "s: " "1256"
## [1] "s: " "1257"
## [1] "s: " "1258"
## [1] "s: " "1259"
## [1] "s: " "1260"
## [1] "s: " "1261"
## [1] "s: " "1262"
## [1] "s: " "1263"
## [1] "s: " "1264"
## [1] "s: " "1265"
## [1] "s: " "1266"
## [1] "s: " "1267"
## [1] "s: " "1268"
## [1] "s: " "1269"
## [1] "s: " "1270"
## [1] "s: " "1271"
## [1] "s: " "1272"
## [1] "s: " "1273"
## [1] "s: " "1274"
## [1] "s: " "1275"
## [1] "s: " "1276"
## [1] "s: " "1277"
## [1] "s: " "1278"
## [1] "s: " "1279"
## [1] "s: " "1280"
## [1] "s: " "1281"
## [1] "s: " "1282"
## [1] "s: " "1283"
## [1] "s: " "1284"
## [1] "s: " "1285"
## [1] "s: " "1286"
## [1] "s: " "1287"
## [1] "s: " "1288"
## [1] "s: " "1289"
## [1] "s: " "1290"
## [1] "s: " "1291"
## [1] "s: " "1292"
## [1] "s: " "1293"
## [1] "s: " "1294"
## [1] "s: " "1295"
## [1] "s: " "1296"
## [1] "s: " "1297"
## [1] "s: " "1298"
## [1] "s: " "1299"
## [1] "s: " "1300"
## [1] "s: " "1301"
## [1] "s: " "1302"
## [1] "s: " "1303"
## [1] "s: " "1304"
## [1] "s: " "1305"
## [1] "s: " "1306"
## [1] "s: " "1307"
## [1] "s: " "1308"
## [1] "s: " "1309"
## [1] "s: " "1310"
## [1] "s: " "1311"
## [1] "s: " "1312"
## [1] "s: " "1313"
## [1] "s: " "1314"
## [1] "s: " "1315"
## [1] "s: " "1316"
## [1] "s: " "1317"
## [1] "s: " "1318"
## [1] "s: " "1319"
## [1] "s: " "1320"
## [1] "s: " "1321"
## [1] "s: " "1322"
## [1] "s: " "1323"
## [1] "s: " "1324"
## [1] "s: " "1325"
## [1] "s: " "1326"
## [1] "s: " "1327"
## [1] "s: " "1328"
## [1] "s: " "1329"
## [1] "s: " "1330"
## [1] "s: " "1331"
## [1] "s: " "1332"
## [1] "s: " "1333"
## [1] "s: " "1334"
## [1] "s: " "1335"
## [1] "s: " "1336"
## [1] "s: " "1337"
## [1] "s: " "1338"
## [1] "s: " "1339"
## [1] "s: " "1340"
## [1] "s: " "1341"
## [1] "s: " "1342"
## [1] "s: " "1343"
## [1] "s: " "1344"
## [1] "s: " "1345"
## [1] "s: " "1346"
## [1] "s: " "1347"
## [1] "s: " "1348"
## [1] "s: " "1349"
## [1] "s: " "1350"
## [1] "s: " "1351"
## [1] "s: " "1352"
## [1] "s: " "1353"
## [1] "s: " "1354"
## [1] "s: " "1355"
## [1] "s: " "1356"
## [1] "s: " "1357"
## [1] "s: " "1358"
## [1] "s: " "1359"
## [1] "s: " "1360"
## [1] "s: " "1361"
## [1] "s: " "1362"
## [1] "s: " "1363"
## [1] "s: " "1364"
## [1] "s: " "1365"
## [1] "s: " "1366"
## [1] "s: " "1367"
## [1] "s: " "1368"
## [1] "s: " "1369"
## [1] "s: " "1370"
## [1] "s: " "1371"
## [1] "s: " "1372"
## [1] "s: " "1373"
## [1] "s: " "1374"
## [1] "s: " "1375"
## [1] "s: " "1376"
## [1] "s: " "1377"
## [1] "s: " "1378"
## [1] "s: " "1379"
## [1] "s: " "1380"
## [1] "s: " "1381"
## [1] "s: " "1382"
## [1] "s: " "1383"
## [1] "s: " "1384"
## [1] "s: " "1385"
## [1] "s: " "1386"
## [1] "s: " "1387"
## [1] "s: " "1388"
## [1] "s: " "1389"
## [1] "s: " "1390"
## [1] "s: " "1391"
## [1] "s: " "1392"
## [1] "s: " "1393"
## [1] "s: " "1394"
## [1] "s: " "1395"
## [1] "s: " "1396"
## [1] "s: " "1397"
## [1] "s: " "1398"
## [1] "s: " "1399"
## [1] "s: " "1400"
## [1] "s: " "1401"
## [1] "s: " "1402"
## [1] "s: " "1403"
## [1] "s: " "1404"
## [1] "s: " "1405"
## [1] "s: " "1406"
## [1] "s: " "1407"
## [1] "s: " "1408"
## [1] "s: " "1409"
## [1] "s: " "1410"
## [1] "s: " "1411"
## [1] "s: " "1412"
## [1] "s: " "1413"
## [1] "s: " "1414"
## [1] "s: " "1415"
## [1] "s: " "1416"
## [1] "s: " "1417"
## [1] "s: " "1418"
## [1] "s: " "1419"
## [1] "s: " "1420"
## [1] "s: " "1421"
## [1] "s: " "1422"
## [1] "s: " "1423"
## [1] "s: " "1424"
## [1] "s: " "1425"
## [1] "s: " "1426"
## [1] "s: " "1427"
## [1] "s: " "1428"
## [1] "s: " "1429"
## [1] "s: " "1430"
## [1] "s: " "1431"
## [1] "s: " "1432"
## [1] "s: " "1433"
## [1] "s: " "1434"
## [1] "s: " "1435"
## [1] "s: " "1436"
## [1] "s: " "1437"
## [1] "s: " "1438"
## [1] "s: " "1439"
## [1] "s: " "1440"
## [1] "s: " "1441"
## [1] "s: " "1442"
## [1] "s: " "1443"
## [1] "s: " "1444"
## [1] "s: " "1445"
## [1] "s: " "1446"
## [1] "s: " "1447"
## [1] "s: " "1448"
## [1] "s: " "1449"
## [1] "s: " "1450"
## [1] "s: " "1451"
## [1] "s: " "1452"
## [1] "s: " "1453"
## [1] "s: " "1454"
## [1] "s: " "1455"
## [1] "s: " "1456"
## [1] "s: " "1457"
## [1] "s: " "1458"
## [1] "s: " "1459"
## [1] "s: " "1460"
## [1] "s: " "1461"
## [1] "s: " "1462"
## [1] "s: " "1463"
## [1] "s: " "1464"
## [1] "s: " "1465"
## [1] "s: " "1466"
## [1] "s: " "1467"
## [1] "s: " "1468"
## [1] "s: " "1469"
## [1] "s: " "1470"
## [1] "s: " "1471"
## [1] "s: " "1472"
## [1] "s: " "1473"
## [1] "s: " "1474"
## [1] "s: " "1475"
## [1] "s: " "1476"
## [1] "s: " "1477"
## [1] "s: " "1478"
## [1] "s: " "1479"
## [1] "s: " "1480"
## [1] "s: " "1481"
## [1] "s: " "1482"
## [1] "s: " "1483"
## [1] "s: " "1484"
## [1] "s: " "1485"
## [1] "s: " "1486"
## [1] "s: " "1487"
## [1] "s: " "1488"
## [1] "s: " "1489"
## [1] "s: " "1490"
## [1] "s: " "1491"
## [1] "s: " "1492"
## [1] "s: " "1493"
## [1] "s: " "1494"
## [1] "s: " "1495"
## [1] "s: " "1496"
## [1] "s: " "1497"
## [1] "s: " "1498"
## [1] "s: " "1499"
## [1] "s: " "1500"
## [1] "s: " "1501"
## [1] "s: " "1502"
## [1] "s: " "1503"
## [1] "s: " "1504"
## [1] "s: " "1505"
## [1] "s: " "1506"
## [1] "s: " "1507"
## [1] "s: " "1508"
## [1] "s: " "1509"
## [1] "s: " "1510"
## [1] "s: " "1511"
## [1] "s: " "1512"
## [1] "s: " "1513"
## [1] "s: " "1514"
## [1] "s: " "1515"
## [1] "s: " "1516"
## [1] "s: " "1517"
## [1] "s: " "1518"
## [1] "s: " "1519"
## [1] "s: " "1520"
## [1] "s: " "1521"
## [1] "s: " "1522"
## [1] "s: " "1523"
## [1] "s: " "1524"
## [1] "s: " "1525"
## [1] "s: " "1526"
## [1] "s: " "1527"
## [1] "s: " "1528"
## [1] "s: " "1529"
## [1] "s: " "1530"
## [1] "s: " "1531"
## [1] "s: " "1532"
## [1] "s: " "1533"
## [1] "s: " "1534"
## [1] "s: " "1535"
## [1] "s: " "1536"
## [1] "s: " "1537"
## [1] "s: " "1538"
## [1] "s: " "1539"
## [1] "s: " "1540"
## [1] "s: " "1541"
## [1] "s: " "1542"
## [1] "s: " "1543"
## [1] "s: " "1544"
## [1] "s: " "1545"
## [1] "s: " "1546"
## [1] "s: " "1547"
## [1] "s: " "1548"
## [1] "s: " "1549"
## [1] "s: " "1550"
## [1] "s: " "1551"
## [1] "s: " "1552"
## [1] "s: " "1553"
## [1] "s: " "1554"
## [1] "s: " "1555"
## [1] "s: " "1556"
## [1] "s: " "1557"
## [1] "s: " "1558"
## [1] "s: " "1559"
## [1] "s: " "1560"
## [1] "s: " "1561"
## [1] "s: " "1562"
## [1] "s: " "1563"
## [1] "s: " "1564"
## [1] "s: " "1565"
## [1] "s: " "1566"
## [1] "s: " "1567"
## [1] "s: " "1568"
## [1] "s: " "1569"
## [1] "s: " "1570"
## [1] "s: " "1571"
## [1] "s: " "1572"
## [1] "s: " "1573"
## [1] "s: " "1574"
## [1] "s: " "1575"
## [1] "s: " "1576"
## [1] "s: " "1577"
## [1] "s: " "1578"
## [1] "s: " "1579"
## [1] "s: " "1580"
## [1] "s: " "1581"
## [1] "s: " "1582"
## [1] "s: " "1583"
## [1] "s: " "1584"
## [1] "s: " "1585"
## [1] "s: " "1586"
## [1] "s: " "1587"
## [1] "s: " "1588"
## [1] "s: " "1589"
## [1] "s: " "1590"
## [1] "s: " "1591"
## [1] "s: " "1592"
## [1] "s: " "1593"
## [1] "s: " "1594"
## [1] "s: " "1595"
## [1] "s: " "1596"
## [1] "s: " "1597"
## [1] "s: " "1598"
## [1] "s: " "1599"
## [1] "s: " "1600"
## [1] "s: " "1601"
## [1] "s: " "1602"
## [1] "s: " "1603"
## [1] "s: " "1604"
## [1] "s: " "1605"
## [1] "s: " "1606"
## [1] "s: " "1607"
## [1] "s: " "1608"
## [1] "s: " "1609"
## [1] "s: " "1610"
## [1] "s: " "1611"
## [1] "s: " "1612"
## [1] "s: " "1613"
## [1] "s: " "1614"
## [1] "s: " "1615"
## [1] "s: " "1616"
## [1] "s: " "1617"
## [1] "s: " "1618"
## [1] "s: " "1619"
## [1] "s: " "1620"
## [1] "s: " "1621"
## [1] "s: " "1622"
## [1] "s: " "1623"
## [1] "s: " "1624"
## [1] "s: " "1625"
## [1] "s: " "1626"
## [1] "s: " "1627"
## [1] "s: " "1628"
## [1] "s: " "1629"
## [1] "s: " "1630"
## [1] "s: " "1631"
## [1] "s: " "1632"
## [1] "s: " "1633"
## [1] "s: " "1634"
## [1] "s: " "1635"
## [1] "s: " "1636"
## [1] "s: " "1637"
## [1] "s: " "1638"
## [1] "s: " "1639"
## [1] "s: " "1640"
## [1] "s: " "1641"
## [1] "s: " "1642"
## [1] "s: " "1643"
## [1] "s: " "1644"
## [1] "s: " "1645"
## [1] "s: " "1646"
## [1] "s: " "1647"
## [1] "s: " "1648"
## [1] "s: " "1649"
## [1] "s: " "1650"
## [1] "s: " "1651"
## [1] "s: " "1652"
## [1] "s: " "1653"
## [1] "s: " "1654"
## [1] "s: " "1655"
## [1] "s: " "1656"
## [1] "s: " "1657"
## [1] "s: " "1658"
## [1] "s: " "1659"
## [1] "s: " "1660"
## [1] "s: " "1661"
## [1] "s: " "1662"
## [1] "s: " "1663"
## [1] "s: " "1664"
## [1] "s: " "1665"
## [1] "s: " "1666"
## [1] "s: " "1667"
## [1] "s: " "1668"
## [1] "s: " "1669"
## [1] "s: " "1670"
## [1] "s: " "1671"
## [1] "s: " "1672"
## [1] "s: " "1673"
## [1] "s: " "1674"
## [1] "s: " "1675"
## [1] "s: " "1676"
## [1] "s: " "1677"
## [1] "s: " "1678"
## [1] "s: " "1679"
## [1] "s: " "1680"
## [1] "s: " "1681"
## [1] "s: " "1682"
## [1] "s: " "1683"
## [1] "s: " "1684"
## [1] "s: " "1685"
## [1] "s: " "1686"
## [1] "s: " "1687"
## [1] "s: " "1688"
## [1] "s: " "1689"
## [1] "s: " "1690"
## [1] "s: " "1691"
## [1] "s: " "1692"
## [1] "s: " "1693"
## [1] "s: " "1694"
## [1] "s: " "1695"
## [1] "s: " "1696"
## [1] "s: " "1697"
## [1] "s: " "1698"
## [1] "s: " "1699"
## [1] "s: " "1700"
## [1] "s: " "1701"
## [1] "s: " "1702"
## [1] "s: " "1703"
## [1] "s: " "1704"
## [1] "s: " "1705"
## [1] "s: " "1706"
## [1] "s: " "1707"
## [1] "s: " "1708"
## [1] "s: " "1709"
## [1] "s: " "1710"
## [1] "s: " "1711"
## [1] "s: " "1712"
## [1] "s: " "1713"
## [1] "s: " "1714"
## [1] "s: " "1715"
## [1] "s: " "1716"
## [1] "s: " "1717"
## [1] "s: " "1718"
## [1] "s: " "1719"
## [1] "s: " "1720"
## [1] "s: " "1721"
## [1] "s: " "1722"
## [1] "s: " "1723"
## [1] "s: " "1724"
## [1] "s: " "1725"
## [1] "s: " "1726"
## [1] "s: " "1727"
## [1] "s: " "1728"
## [1] "s: " "1729"
## [1] "s: " "1730"
## [1] "s: " "1731"
## [1] "s: " "1732"
## [1] "s: " "1733"
## [1] "s: " "1734"
## [1] "s: " "1735"
## [1] "s: " "1736"
## [1] "s: " "1737"
## [1] "s: " "1738"
## [1] "s: " "1739"
## [1] "s: " "1740"
## [1] "s: " "1741"
## [1] "s: " "1742"
## [1] "s: " "1743"
## [1] "s: " "1744"
## [1] "s: " "1745"
## [1] "s: " "1746"
## [1] "s: " "1747"
## [1] "s: " "1748"
## [1] "s: " "1749"
## [1] "s: " "1750"
## [1] "s: " "1751"
## [1] "s: " "1752"
## [1] "s: " "1753"
## [1] "s: " "1754"
## [1] "s: " "1755"
## [1] "s: " "1756"
## [1] "s: " "1757"
## [1] "s: " "1758"
## [1] "s: " "1759"
## [1] "s: " "1760"
## [1] "s: " "1761"
## [1] "s: " "1762"
## [1] "s: " "1763"
## [1] "s: " "1764"
## [1] "s: " "1765"
## [1] "s: " "1766"
## [1] "s: " "1767"
## [1] "s: " "1768"
## [1] "s: " "1769"
## [1] "s: " "1770"
## [1] "s: " "1771"
## [1] "s: " "1772"
## [1] "s: " "1773"
## [1] "s: " "1774"
## [1] "s: " "1775"
## [1] "s: " "1776"
## [1] "s: " "1777"
## [1] "s: " "1778"
## [1] "s: " "1779"
## [1] "s: " "1780"
## [1] "s: " "1781"
## [1] "s: " "1782"
## [1] "s: " "1783"
## [1] "s: " "1784"
## [1] "s: " "1785"
## [1] "s: " "1786"
## [1] "s: " "1787"
## [1] "s: " "1788"
## [1] "s: " "1789"
## [1] "s: " "1790"
## [1] "s: " "1791"
## [1] "s: " "1792"
## [1] "s: " "1793"
## [1] "s: " "1794"
## [1] "s: " "1795"
## [1] "s: " "1796"
## [1] "s: " "1797"
## [1] "s: " "1798"
## [1] "s: " "1799"
## [1] "s: " "1800"
## [1] "s: " "1801"
## [1] "s: " "1802"
## [1] "s: " "1803"
## [1] "s: " "1804"
## [1] "s: " "1805"
## [1] "s: " "1806"
## [1] "s: " "1807"
## [1] "s: " "1808"
## [1] "s: " "1809"
## [1] "s: " "1810"
## [1] "s: " "1811"
## [1] "s: " "1812"
## [1] "s: " "1813"
## [1] "s: " "1814"
## [1] "s: " "1815"
## [1] "s: " "1816"
## [1] "s: " "1817"
## [1] "s: " "1818"
## [1] "s: " "1819"
## [1] "s: " "1820"
## [1] "s: " "1821"
## [1] "s: " "1822"
## [1] "s: " "1823"
## [1] "s: " "1824"
## [1] "s: " "1825"
## [1] "s: " "1826"
## [1] "s: " "1827"
## [1] "s: " "1828"
## [1] "s: " "1829"
## [1] "s: " "1830"
## [1] "s: " "1831"
## [1] "s: " "1832"
## [1] "s: " "1833"
## [1] "s: " "1834"
## [1] "s: " "1835"
## [1] "s: " "1836"
## [1] "s: " "1837"
## [1] "s: " "1838"
## [1] "s: " "1839"
## [1] "s: " "1840"
## [1] "s: " "1841"
## [1] "s: " "1842"
## [1] "s: " "1843"
## [1] "s: " "1844"
## [1] "s: " "1845"
## [1] "s: " "1846"
## [1] "s: " "1847"
## [1] "s: " "1848"
## [1] "s: " "1849"
## [1] "s: " "1850"
## [1] "s: " "1851"
## [1] "s: " "1852"
## [1] "s: " "1853"
## [1] "s: " "1854"
## [1] "s: " "1855"
## [1] "s: " "1856"
## [1] "s: " "1857"
## [1] "s: " "1858"
## [1] "s: " "1859"
## [1] "s: " "1860"
## [1] "s: " "1861"
## [1] "s: " "1862"
## [1] "s: " "1863"
## [1] "s: " "1864"
## [1] "s: " "1865"
## [1] "s: " "1866"
## [1] "s: " "1867"
## [1] "s: " "1868"
## [1] "s: " "1869"
## [1] "s: " "1870"
## [1] "s: " "1871"
## [1] "s: " "1872"
## [1] "s: " "1873"
## [1] "s: " "1874"
## [1] "s: " "1875"
## [1] "s: " "1876"
## [1] "s: " "1877"
## [1] "s: " "1878"
## [1] "s: " "1879"
## [1] "s: " "1880"
## [1] "s: " "1881"
## [1] "s: " "1882"
## [1] "s: " "1883"
## [1] "s: " "1884"
## [1] "s: " "1885"
## [1] "s: " "1886"
## [1] "s: " "1887"
## [1] "s: " "1888"
## [1] "s: " "1889"
## [1] "s: " "1890"
## [1] "s: " "1891"
## [1] "s: " "1892"
## [1] "s: " "1893"
## [1] "s: " "1894"
## [1] "s: " "1895"
## [1] "s: " "1896"
## [1] "s: " "1897"
## [1] "s: " "1898"
## [1] "s: " "1899"
## [1] "s: " "1900"
## [1] "s: " "1901"
## [1] "s: " "1902"
## [1] "s: " "1903"
## [1] "s: " "1904"
## [1] "s: " "1905"
## [1] "s: " "1906"
## [1] "s: " "1907"
## [1] "s: " "1908"
## [1] "s: " "1909"
## [1] "s: " "1910"
## [1] "s: " "1911"
## [1] "s: " "1912"
## [1] "s: " "1913"
## [1] "s: " "1914"
## [1] "s: " "1915"
## [1] "s: " "1916"
## [1] "s: " "1917"
## [1] "s: " "1918"
## [1] "s: " "1919"
## [1] "s: " "1920"
## [1] "s: " "1921"
## [1] "s: " "1922"
## [1] "s: " "1923"
## [1] "s: " "1924"
## [1] "s: " "1925"
## [1] "s: " "1926"
## [1] "s: " "1927"
## [1] "s: " "1928"
## [1] "s: " "1929"
## [1] "s: " "1930"
## [1] "s: " "1931"
## [1] "s: " "1932"
## [1] "s: " "1933"
## [1] "s: " "1934"
## [1] "s: " "1935"
## [1] "s: " "1936"
## [1] "s: " "1937"
## [1] "s: " "1938"
## [1] "s: " "1939"
## [1] "s: " "1940"
## [1] "s: " "1941"
## [1] "s: " "1942"
## [1] "s: " "1943"
## [1] "s: " "1944"
## [1] "s: " "1945"
## [1] "s: " "1946"
## [1] "s: " "1947"
## [1] "s: " "1948"
## [1] "s: " "1949"
## [1] "s: " "1950"
## [1] "s: " "1951"
## [1] "s: " "1952"
## [1] "s: " "1953"
## [1] "s: " "1954"
## [1] "s: " "1955"
## [1] "s: " "1956"
## [1] "s: " "1957"
## [1] "s: " "1958"
## [1] "s: " "1959"
## [1] "s: " "1960"
## [1] "s: " "1961"
## [1] "s: " "1962"
## [1] "s: " "1963"
## [1] "s: " "1964"
## [1] "s: " "1965"
## [1] "s: " "1966"
## [1] "s: " "1967"
## [1] "s: " "1968"
## [1] "s: " "1969"
## [1] "s: " "1970"
## [1] "s: " "1971"
## [1] "s: " "1972"
## [1] "s: " "1973"
## [1] "s: " "1974"
## [1] "s: " "1975"
## [1] "s: " "1976"
## [1] "s: " "1977"
## [1] "s: " "1978"
## [1] "s: " "1979"
## [1] "s: " "1980"
## [1] "s: " "1981"
## [1] "s: " "1982"
## [1] "s: " "1983"
## [1] "s: " "1984"
## [1] "s: " "1985"
## [1] "s: " "1986"
## [1] "s: " "1987"
## [1] "s: " "1988"
## [1] "s: " "1989"
## [1] "s: " "1990"
## [1] "s: " "1991"
## [1] "s: " "1992"
## [1] "s: " "1993"
## [1] "s: " "1994"
## [1] "s: " "1995"
## [1] "s: " "1996"
## [1] "s: " "1997"
## [1] "s: " "1998"
## [1] "s: " "1999"
## [1] "s: " "2000"
## [1] "s: " "2001"
## [1] "s: " "2002"
## [1] "s: " "2003"
## [1] "s: " "2004"
## [1] "s: " "2005"
## [1] "s: " "2006"
## [1] "s: " "2007"
## [1] "s: " "2008"
## [1] "s: " "2009"
## [1] "s: " "2010"
## [1] "s: " "2011"
## [1] "s: " "2012"
## [1] "s: " "2013"
## [1] "s: " "2014"
## [1] "s: " "2015"
## [1] "s: " "2016"
## [1] "s: " "2017"
## [1] "s: " "2018"
## [1] "s: " "2019"
## [1] "s: " "2020"
## [1] "s: " "2021"
## [1] "s: " "2022"
## [1] "s: " "2023"
## [1] "s: " "2024"
## [1] "s: " "2025"
## [1] "s: " "2026"
## [1] "s: " "2027"
## [1] "s: " "2028"
## [1] "s: " "2029"
## [1] "s: " "2030"
## [1] "s: " "2031"
## [1] "s: " "2032"
## [1] "s: " "2033"
## [1] "s: " "2034"
## [1] "s: " "2035"
## [1] "s: " "2036"
## [1] "s: " "2037"
## [1] "s: " "2038"
## [1] "s: " "2039"
## [1] "s: " "2040"
## [1] "s: " "2041"
## [1] "s: " "2042"
## [1] "s: " "2043"
## [1] "s: " "2044"
## [1] "s: " "2045"
## [1] "s: " "2046"
## [1] "s: " "2047"
## [1] "s: " "2048"
## [1] "s: " "2049"
## [1] "s: " "2050"
## [1] "s: " "2051"
## [1] "s: " "2052"
## [1] "s: " "2053"
## [1] "s: " "2054"
## [1] "s: " "2055"
## [1] "s: " "2056"
## [1] "s: " "2057"
## [1] "s: " "2058"
## [1] "s: " "2059"
## [1] "s: " "2060"
## [1] "s: " "2061"
## [1] "s: " "2062"
## [1] "s: " "2063"
## [1] "s: " "2064"
## [1] "s: " "2065"
## [1] "s: " "2066"
## [1] "s: " "2067"
## [1] "s: " "2068"
## [1] "s: " "2069"
## [1] "s: " "2070"
## [1] "s: " "2071"
## [1] "s: " "2072"
## [1] "s: " "2073"
## [1] "s: " "2074"
## [1] "s: " "2075"
## [1] "s: " "2076"
## [1] "s: " "2077"
## [1] "s: " "2078"
## [1] "s: " "2079"
## [1] "s: " "2080"
## [1] "s: " "2081"
## [1] "s: " "2082"
## [1] "s: " "2083"
## [1] "s: " "2084"
## [1] "s: " "2085"
## [1] "s: " "2086"
## [1] "s: " "2087"
## [1] "s: " "2088"
## [1] "s: " "2089"
## [1] "s: " "2090"
## [1] "s: " "2091"
## [1] "s: " "2092"
## [1] "s: " "2093"
## [1] "s: " "2094"
## [1] "s: " "2095"
## [1] "s: " "2096"
## [1] "s: " "2097"
## [1] "s: " "2098"
## [1] "s: " "2099"
## [1] "s: " "2100"
## [1] "s: " "2101"
## [1] "s: " "2102"
## [1] "s: " "2103"
## [1] "s: " "2104"
## [1] "s: " "2105"
## [1] "s: " "2106"
## [1] "s: " "2107"
## [1] "s: " "2108"
## [1] "s: " "2109"
## [1] "s: " "2110"
## [1] "s: " "2111"
## [1] "s: " "2112"
## [1] "s: " "2113"
## [1] "s: " "2114"
## [1] "s: " "2115"
## [1] "s: " "2116"
## [1] "s: " "2117"
## [1] "s: " "2118"
## [1] "s: " "2119"
## [1] "s: " "2120"
## [1] "s: " "2121"
## [1] "s: " "2122"
## [1] "s: " "2123"
## [1] "s: " "2124"
## [1] "s: " "2125"
## [1] "s: " "2126"
## [1] "s: " "2127"
## [1] "s: " "2128"
## [1] "s: " "2129"
## [1] "s: " "2130"
## [1] "s: " "2131"
## [1] "s: " "2132"
## [1] "s: " "2133"
## [1] "s: " "2134"
## [1] "s: " "2135"
## [1] "s: " "2136"
## [1] "s: " "2137"
## [1] "s: " "2138"
## [1] "s: " "2139"
## [1] "s: " "2140"
## [1] "s: " "2141"
## [1] "s: " "2142"
## [1] "s: " "2143"
## [1] "s: " "2144"
## [1] "s: " "2145"
## [1] "s: " "2146"
## [1] "s: " "2147"
## [1] "s: " "2148"
## [1] "s: " "2149"
## [1] "s: " "2150"
## [1] "s: " "2151"
## [1] "s: " "2152"
## [1] "s: " "2153"
## [1] "s: " "2154"
## [1] "s: " "2155"
## [1] "s: " "2156"
## [1] "s: " "2157"
## [1] "s: " "2158"
## [1] "s: " "2159"
## [1] "s: " "2160"
## [1] "s: " "2161"
## [1] "s: " "2162"
## [1] "s: " "2163"
## [1] "s: " "2164"
## [1] "s: " "2165"
## [1] "s: " "2166"
## [1] "s: " "2167"
## [1] "s: " "2168"
## [1] "s: " "2169"
## [1] "s: " "2170"
## [1] "s: " "2171"
## [1] "s: " "2172"
## [1] "s: " "2173"
## [1] "s: " "2174"
## [1] "s: " "2175"
## [1] "s: " "2176"
## [1] "s: " "2177"
## [1] "s: " "2178"
## [1] "s: " "2179"
## [1] "s: " "2180"
## [1] "s: " "2181"
## [1] "s: " "2182"
## [1] "s: " "2183"
## [1] "s: " "2184"
## [1] "s: " "2185"
## [1] "s: " "2186"
## [1] "s: " "2187"
## [1] "s: " "2188"
## [1] "s: " "2189"
## [1] "s: " "2190"
## [1] "s: " "2191"
## [1] "s: " "2192"
## [1] "s: " "2193"
## [1] "s: " "2194"
## [1] "s: " "2195"
## [1] "s: " "2196"
## [1] "s: " "2197"
## [1] "s: " "2198"
## [1] "s: " "2199"
## [1] "s: " "2200"
## [1] "s: " "2201"
## [1] "s: " "2202"
## [1] "s: " "2203"
## [1] "s: " "2204"
## [1] "s: " "2205"
## [1] "s: " "2206"
## [1] "s: " "2207"
## [1] "s: " "2208"
## [1] "s: " "2209"
## [1] "s: " "2210"
## [1] "s: " "2211"
## [1] "s: " "2212"
## [1] "s: " "2213"
## [1] "s: " "2214"
## [1] "s: " "2215"
## [1] "s: " "2216"
## [1] "s: " "2217"
## [1] "s: " "2218"
## [1] "s: " "2219"
## [1] "s: " "2220"
## [1] "s: " "2221"
## [1] "s: " "2222"
## [1] "s: " "2223"
## [1] "s: " "2224"
## [1] "s: " "2225"
## [1] "s: " "2226"
## [1] "s: " "2227"
## [1] "s: " "2228"
## [1] "s: " "2229"
## [1] "s: " "2230"
## [1] "s: " "2231"
## [1] "s: " "2232"
## [1] "s: " "2233"
## [1] "s: " "2234"
## [1] "s: " "2235"
## [1] "s: " "2236"
## [1] "s: " "2237"
## [1] "s: " "2238"
## [1] "s: " "2239"
## [1] "s: " "2240"
## [1] "s: " "2241"
## [1] "s: " "2242"
## [1] "s: " "2243"
## [1] "s: " "2244"
## [1] "s: " "2245"
## [1] "s: " "2246"
## [1] "s: " "2247"
## [1] "s: " "2248"
## [1] "s: " "2249"
## [1] "s: " "2250"
## [1] "s: " "2251"
## [1] "s: " "2252"
## [1] "s: " "2253"
## [1] "s: " "2254"
## [1] "s: " "2255"
## [1] "s: " "2256"
## [1] "s: " "2257"
## [1] "s: " "2258"
## [1] "s: " "2259"
## [1] "s: " "2260"
## [1] "s: " "2261"
## [1] "s: " "2262"
## [1] "s: " "2263"
## [1] "s: " "2264"
## [1] "s: " "2265"
## [1] "s: " "2266"
## [1] "s: " "2267"
## [1] "s: " "2268"
## [1] "s: " "2269"
## [1] "s: " "2270"
## [1] "s: " "2271"
## [1] "s: " "2272"
## [1] "s: " "2273"
## [1] "s: " "2274"
## [1] "s: " "2275"
## [1] "s: " "2276"
## [1] "s: " "2277"
## [1] "s: " "2278"
## [1] "s: " "2279"
## [1] "s: " "2280"
## [1] "s: " "2281"
## [1] "s: " "2282"
## [1] "s: " "2283"
## [1] "s: " "2284"
## [1] "s: " "2285"
## [1] "s: " "2286"
## [1] "s: " "2287"
## [1] "s: " "2288"
## [1] "s: " "2289"
## [1] "s: " "2290"
## [1] "s: " "2291"
## [1] "s: " "2292"
## [1] "s: " "2293"
## [1] "s: " "2294"
## [1] "s: " "2295"
## [1] "s: " "2296"
## [1] "s: " "2297"
## [1] "s: " "2298"
## [1] "s: " "2299"
## [1] "s: " "2300"
## [1] "s: " "2301"
## [1] "s: " "2302"
## [1] "s: " "2303"
## [1] "s: " "2304"
## [1] "s: " "2305"
## [1] "s: " "2306"
## [1] "s: " "2307"
## [1] "s: " "2308"
## [1] "s: " "2309"
## [1] "s: " "2310"
## [1] "s: " "2311"
## [1] "s: " "2312"
## [1] "s: " "2313"
## [1] "s: " "2314"
## [1] "s: " "2315"
## [1] "s: " "2316"
## [1] "s: " "2317"
## [1] "s: " "2318"
## [1] "s: " "2319"
## [1] "s: " "2320"
## [1] "s: " "2321"
## [1] "s: " "2322"
## [1] "s: " "2323"
## [1] "s: " "2324"
## [1] "s: " "2325"
## [1] "s: " "2326"
## [1] "s: " "2327"
## [1] "s: " "2328"
## [1] "s: " "2329"
## [1] "s: " "2330"
## [1] "s: " "2331"
## [1] "s: " "2332"
## [1] "s: " "2333"
## [1] "s: " "2334"
## [1] "s: " "2335"
## [1] "s: " "2336"
## [1] "s: " "2337"
## [1] "s: " "2338"
## [1] "s: " "2339"
## [1] "s: " "2340"
## [1] "s: " "2341"
## [1] "s: " "2342"
## [1] "s: " "2343"
## [1] "s: " "2344"
## [1] "s: " "2345"
## [1] "s: " "2346"
## [1] "s: " "2347"
## [1] "s: " "2348"
## [1] "s: " "2349"
## [1] "s: " "2350"
## [1] "s: " "2351"
## [1] "s: " "2352"
## [1] "s: " "2353"
## [1] "s: " "2354"
## [1] "s: " "2355"
## [1] "s: " "2356"
## [1] "s: " "2357"
## [1] "s: " "2358"
## [1] "s: " "2359"
## [1] "s: " "2360"
## [1] "s: " "2361"
## [1] "s: " "2362"
## [1] "s: " "2363"
## [1] "s: " "2364"
## [1] "s: " "2365"
## [1] "s: " "2366"
## [1] "s: " "2367"
## [1] "s: " "2368"
## [1] "s: " "2369"
## [1] "s: " "2370"
## [1] "s: " "2371"
## [1] "s: " "2372"
## [1] "s: " "2373"
## [1] "s: " "2374"
## [1] "s: " "2375"
## [1] "s: " "2376"
## [1] "s: " "2377"
## [1] "s: " "2378"
## [1] "s: " "2379"
## [1] "s: " "2380"
## [1] "s: " "2381"
## [1] "s: " "2382"
## [1] "s: " "2383"
## [1] "s: " "2384"
## [1] "s: " "2385"
## [1] "s: " "2386"
## [1] "s: " "2387"
## [1] "s: " "2388"
## [1] "s: " "2389"
## [1] "s: " "2390"
## [1] "s: " "2391"
## [1] "s: " "2392"
## [1] "s: " "2393"
## [1] "s: " "2394"
## [1] "s: " "2395"
## [1] "s: " "2396"
## [1] "s: " "2397"
## [1] "s: " "2398"
## [1] "s: " "2399"
## [1] "s: " "2400"
## [1] "s: " "2401"
## [1] "s: " "2402"
## [1] "s: " "2403"
## [1] "s: " "2404"
## [1] "s: " "2405"
## [1] "s: " "2406"
## [1] "s: " "2407"
## [1] "s: " "2408"
## [1] "s: " "2409"
## [1] "s: " "2410"
## [1] "s: " "2411"
## [1] "s: " "2412"
## [1] "s: " "2413"
## [1] "s: " "2414"
## [1] "s: " "2415"
## [1] "s: " "2416"
## [1] "s: " "2417"
## [1] "s: " "2418"
## [1] "s: " "2419"
## [1] "s: " "2420"
## [1] "s: " "2421"
## [1] "s: " "2422"
## [1] "s: " "2423"
## [1] "s: " "2424"
## [1] "s: " "2425"
## [1] "s: " "2426"
## [1] "s: " "2427"
## [1] "s: " "2428"
## [1] "s: " "2429"
## [1] "s: " "2430"
## [1] "s: " "2431"
## [1] "s: " "2432"
## [1] "s: " "2433"
## [1] "s: " "2434"
## [1] "s: " "2435"
## [1] "s: " "2436"
## [1] "s: " "2437"
## [1] "s: " "2438"
## [1] "s: " "2439"
## [1] "s: " "2440"
## [1] "s: " "2441"
## [1] "s: " "2442"
## [1] "s: " "2443"
## [1] "s: " "2444"
## [1] "s: " "2445"
## [1] "s: " "2446"
## [1] "s: " "2447"
## [1] "s: " "2448"
## [1] "s: " "2449"
## [1] "s: " "2450"
## [1] "s: " "2451"
## [1] "s: " "2452"
## [1] "s: " "2453"
## [1] "s: " "2454"
## [1] "s: " "2455"
## [1] "s: " "2456"
## [1] "s: " "2457"
## [1] "s: " "2458"
## [1] "s: " "2459"
## [1] "s: " "2460"
## [1] "s: " "2461"
## [1] "s: " "2462"
## [1] "s: " "2463"
## [1] "s: " "2464"
## [1] "s: " "2465"
## [1] "s: " "2466"
## [1] "s: " "2467"
## [1] "s: " "2468"
## [1] "s: " "2469"
## [1] "s: " "2470"
## [1] "s: " "2471"
## [1] "s: " "2472"
## [1] "s: " "2473"
## [1] "s: " "2474"
## [1] "s: " "2475"
## [1] "s: " "2476"
## [1] "s: " "2477"
## [1] "s: " "2478"
## [1] "s: " "2479"
## [1] "s: " "2480"
## [1] "s: " "2481"
## [1] "s: " "2482"
## [1] "s: " "2483"
## [1] "s: " "2484"
## [1] "s: " "2485"
## [1] "s: " "2486"
## [1] "s: " "2487"
## [1] "s: " "2488"
## [1] "s: " "2489"
## [1] "s: " "2490"
## [1] "s: " "2491"
## [1] "s: " "2492"
## [1] "s: " "2493"
## [1] "s: " "2494"
## [1] "s: " "2495"
## [1] "s: " "2496"
## [1] "s: " "2497"
## [1] "s: " "2498"
## [1] "s: " "2499"
## [1] "s: " "2500"
## [1] "s: " "2501"
## [1] "s: " "2502"
## [1] "s: " "2503"
## [1] "s: " "2504"
## [1] "s: " "2505"
## [1] "s: " "2506"
## [1] "s: " "2507"
## [1] "s: " "2508"
## [1] "s: " "2509"
## [1] "s: " "2510"
## [1] "s: " "2511"
## [1] "s: " "2512"
## [1] "s: " "2513"
## [1] "s: " "2514"
## [1] "s: " "2515"
## [1] "s: " "2516"
## [1] "s: " "2517"
## [1] "s: " "2518"
## [1] "s: " "2519"
## [1] "s: " "2520"
## [1] "s: " "2521"
## [1] "s: " "2522"
## [1] "s: " "2523"
## [1] "s: " "2524"
## [1] "s: " "2525"
## [1] "s: " "2526"
## [1] "s: " "2527"
## [1] "s: " "2528"
## [1] "s: " "2529"
## [1] "s: " "2530"
## [1] "s: " "2531"
## [1] "s: " "2532"
## [1] "s: " "2533"
## [1] "s: " "2534"
## [1] "s: " "2535"
## [1] "s: " "2536"
## [1] "s: " "2537"
## [1] "s: " "2538"
## [1] "s: " "2539"
## [1] "s: " "2540"
## [1] "s: " "2541"
## [1] "s: " "2542"
## [1] "s: " "2543"
## [1] "s: " "2544"
## [1] "s: " "2545"
## [1] "s: " "2546"
## [1] "s: " "2547"
## [1] "s: " "2548"
## [1] "s: " "2549"
## [1] "s: " "2550"
## [1] "s: " "2551"
## [1] "s: " "2552"
## [1] "s: " "2553"
## [1] "s: " "2554"
## [1] "s: " "2555"
## [1] "s: " "2556"
## [1] "s: " "2557"
## [1] "s: " "2558"
## [1] "s: " "2559"
## [1] "s: " "2560"
## [1] "s: " "2561"
## [1] "s: " "2562"
## [1] "s: " "2563"
## [1] "s: " "2564"
## [1] "s: " "2565"
## [1] "s: " "2566"
## [1] "s: " "2567"
## [1] "s: " "2568"
## [1] "s: " "2569"
## [1] "s: " "2570"
## [1] "s: " "2571"
## [1] "s: " "2572"
## [1] "s: " "2573"
## [1] "s: " "2574"
## [1] "s: " "2575"
## [1] "s: " "2576"
## [1] "s: " "2577"
## [1] "s: " "2578"
## [1] "s: " "2579"
## [1] "s: " "2580"
## [1] "s: " "2581"
## [1] "s: " "2582"
## [1] "s: " "2583"
## [1] "s: " "2584"
## [1] "s: " "2585"
## [1] "s: " "2586"
## [1] "s: " "2587"
## [1] "s: " "2588"
## [1] "s: " "2589"
## [1] "s: " "2590"
## [1] "s: " "2591"
## [1] "s: " "2592"
## [1] "s: " "2593"
## [1] "s: " "2594"
## [1] "s: " "2595"
## [1] "s: " "2596"
## [1] "s: " "2597"
## [1] "s: " "2598"
## [1] "s: " "2599"
## [1] "s: " "2600"
## [1] "s: " "2601"
## [1] "s: " "2602"
## [1] "s: " "2603"
## [1] "s: " "2604"
## [1] "s: " "2605"
## [1] "s: " "2606"
## [1] "s: " "2607"
## [1] "s: " "2608"
## [1] "s: " "2609"
## [1] "s: " "2610"
## [1] "s: " "2611"
## [1] "s: " "2612"
## [1] "s: " "2613"
## [1] "s: " "2614"
## [1] "s: " "2615"
## [1] "s: " "2616"
## [1] "s: " "2617"
## [1] "s: " "2618"
## [1] "s: " "2619"
## [1] "s: " "2620"
## [1] "s: " "2621"
## [1] "s: " "2622"
## [1] "s: " "2623"
## [1] "s: " "2624"
## [1] "s: " "2625"
## [1] "s: " "2626"
## [1] "s: " "2627"
## [1] "s: " "2628"
## [1] "s: " "2629"
## [1] "s: " "2630"
## [1] "s: " "2631"
## [1] "s: " "2632"
## [1] "s: " "2633"
## [1] "s: " "2634"
## [1] "s: " "2635"
## [1] "s: " "2636"
## [1] "s: " "2637"
## [1] "s: " "2638"
## [1] "s: " "2639"
## [1] "s: " "2640"
## [1] "s: " "2641"
## [1] "s: " "2642"
## [1] "s: " "2643"
## [1] "s: " "2644"
## [1] "s: " "2645"
## [1] "s: " "2646"
## [1] "s: " "2647"
## [1] "s: " "2648"
## [1] "s: " "2649"
## [1] "s: " "2650"
## [1] "s: " "2651"
## [1] "s: " "2652"
## [1] "s: " "2653"
## [1] "s: " "2654"
## [1] "s: " "2655"
## [1] "s: " "2656"
## [1] "s: " "2657"
## [1] "s: " "2658"
## [1] "s: " "2659"
## [1] "s: " "2660"
## [1] "s: " "2661"
## [1] "s: " "2662"
## [1] "s: " "2663"
## [1] "s: " "2664"
## [1] "s: " "2665"
## [1] "s: " "2666"
## [1] "s: " "2667"
## [1] "s: " "2668"
## [1] "s: " "2669"
## [1] "s: " "2670"
## [1] "s: " "2671"
## [1] "s: " "2672"
## [1] "s: " "2673"
## [1] "s: " "2674"
## [1] "s: " "2675"
## [1] "s: " "2676"
## [1] "s: " "2677"
## [1] "s: " "2678"
## [1] "s: " "2679"
## [1] "s: " "2680"
## [1] "s: " "2681"
## [1] "s: " "2682"
## [1] "s: " "2683"
## [1] "s: " "2684"
## [1] "s: " "2685"
## [1] "s: " "2686"
## [1] "s: " "2687"
## [1] "s: " "2688"
## [1] "s: " "2689"
## [1] "s: " "2690"
## [1] "s: " "2691"
## [1] "s: " "2692"
## [1] "s: " "2693"
## [1] "s: " "2694"
## [1] "s: " "2695"
## [1] "s: " "2696"
## [1] "s: " "2697"
## [1] "s: " "2698"
## [1] "s: " "2699"
## [1] "s: " "2700"
## [1] "s: " "2701"
## [1] "s: " "2702"
## [1] "s: " "2703"
## [1] "s: " "2704"
## [1] "s: " "2705"
## [1] "s: " "2706"
## [1] "s: " "2707"
## [1] "s: " "2708"
## [1] "s: " "2709"
## [1] "s: " "2710"
## [1] "s: " "2711"
## [1] "s: " "2712"
## [1] "s: " "2713"
## [1] "s: " "2714"
## [1] "s: " "2715"
## [1] "s: " "2716"
## [1] "s: " "2717"
## [1] "s: " "2718"
## [1] "s: " "2719"
## [1] "s: " "2720"
## [1] "s: " "2721"
## [1] "s: " "2722"
## [1] "s: " "2723"
## [1] "s: " "2724"
## [1] "s: " "2725"
## [1] "s: " "2726"
## [1] "s: " "2727"
## [1] "s: " "2728"
## [1] "s: " "2729"
## [1] "s: " "2730"
## [1] "s: " "2731"
## [1] "s: " "2732"
## [1] "s: " "2733"
## [1] "s: " "2734"
## [1] "s: " "2735"
## [1] "s: " "2736"
## [1] "s: " "2737"
## [1] "s: " "2738"
## [1] "s: " "2739"
## [1] "s: " "2740"
## [1] "s: " "2741"
## [1] "s: " "2742"
## [1] "s: " "2743"
## [1] "s: " "2744"
## [1] "s: " "2745"
## [1] "s: " "2746"
## [1] "s: " "2747"
## [1] "s: " "2748"
## [1] "s: " "2749"
## [1] "s: " "2750"
## [1] "s: " "2751"
## [1] "s: " "2752"
## [1] "s: " "2753"
## [1] "s: " "2754"
## [1] "s: " "2755"
## [1] "s: " "2756"
## [1] "s: " "2757"
## [1] "s: " "2758"
## [1] "s: " "2759"
## [1] "s: " "2760"
## [1] "s: " "2761"
## [1] "s: " "2762"
## [1] "s: " "2763"
## [1] "s: " "2764"
## [1] "s: " "2765"
## [1] "s: " "2766"
## [1] "s: " "2767"
## [1] "s: " "2768"
## [1] "s: " "2769"
## [1] "s: " "2770"
## [1] "s: " "2771"
## [1] "s: " "2772"
## [1] "s: " "2773"
## [1] "s: " "2774"
## [1] "s: " "2775"
## [1] "s: " "2776"
## [1] "s: " "2777"
## [1] "s: " "2778"
## [1] "s: " "2779"
## [1] "s: " "2780"
## [1] "s: " "2781"
## [1] "s: " "2782"
## [1] "s: " "2783"
## [1] "s: " "2784"
## [1] "s: " "2785"
## [1] "s: " "2786"
## [1] "s: " "2787"
## [1] "s: " "2788"
## [1] "s: " "2789"
## [1] "s: " "2790"
## [1] "s: " "2791"
## [1] "s: " "2792"
## [1] "s: " "2793"
## [1] "s: " "2794"
## [1] "s: " "2795"
## [1] "s: " "2796"
## [1] "s: " "2797"
## [1] "s: " "2798"
## [1] "s: " "2799"
## [1] "s: " "2800"
## [1] "s: " "2801"
## [1] "s: " "2802"
## [1] "s: " "2803"
## [1] "s: " "2804"
## [1] "s: " "2805"
## [1] "s: " "2806"
## [1] "s: " "2807"
## [1] "s: " "2808"
## [1] "s: " "2809"
## [1] "s: " "2810"
## [1] "s: " "2811"
## [1] "s: " "2812"
## [1] "s: " "2813"
## [1] "s: " "2814"
## [1] "s: " "2815"
## [1] "s: " "2816"
## [1] "s: " "2817"
## [1] "s: " "2818"
## [1] "s: " "2819"
## [1] "s: " "2820"
## [1] "s: " "2821"
## [1] "s: " "2822"
## [1] "s: " "2823"
## [1] "s: " "2824"
## [1] "s: " "2825"
## [1] "s: " "2826"
## [1] "s: " "2827"
## [1] "s: " "2828"
## [1] "s: " "2829"
## [1] "s: " "2830"
## [1] "s: " "2831"
## [1] "s: " "2832"
## [1] "s: " "2833"
## [1] "s: " "2834"
## [1] "s: " "2835"
## [1] "s: " "2836"
## [1] "s: " "2837"
## [1] "s: " "2838"
## [1] "s: " "2839"
## [1] "s: " "2840"
## [1] "s: " "2841"
## [1] "s: " "2842"
## [1] "s: " "2843"
## [1] "s: " "2844"
## [1] "s: " "2845"
## [1] "s: " "2846"
## [1] "s: " "2847"
## [1] "s: " "2848"
## [1] "s: " "2849"
## [1] "s: " "2850"
## [1] "s: " "2851"
## [1] "s: " "2852"
## [1] "s: " "2853"
## [1] "s: " "2854"
## [1] "s: " "2855"
## [1] "s: " "2856"
## [1] "s: " "2857"
## [1] "s: " "2858"
## [1] "s: " "2859"
## [1] "s: " "2860"
## [1] "s: " "2861"
## [1] "s: " "2862"
## [1] "s: " "2863"
## [1] "s: " "2864"
## [1] "s: " "2865"
## [1] "s: " "2866"
## [1] "s: " "2867"
## [1] "s: " "2868"
## [1] "s: " "2869"
## [1] "s: " "2870"
## [1] "s: " "2871"
## [1] "s: " "2872"
## [1] "s: " "2873"
## [1] "s: " "2874"
## [1] "s: " "2875"
## [1] "s: " "2876"
## [1] "s: " "2877"
## [1] "s: " "2878"
## [1] "s: " "2879"
## [1] "s: " "2880"
## [1] "s: " "2881"
## [1] "s: " "2882"
## [1] "s: " "2883"
## [1] "s: " "2884"
## [1] "s: " "2885"
## [1] "s: " "2886"
## [1] "s: " "2887"
## [1] "s: " "2888"
## [1] "s: " "2889"
## [1] "s: " "2890"
## [1] "s: " "2891"
## [1] "s: " "2892"
## [1] "s: " "2893"
## [1] "s: " "2894"
## [1] "s: " "2895"
## [1] "s: " "2896"
## [1] "s: " "2897"
## [1] "s: " "2898"
## [1] "s: " "2899"
## [1] "s: " "2900"
## [1] "s: " "2901"
## [1] "s: " "2902"
## [1] "s: " "2903"
## [1] "s: " "2904"
## [1] "s: " "2905"
## [1] "s: " "2906"
## [1] "s: " "2907"
## [1] "s: " "2908"
## [1] "s: " "2909"
## [1] "s: " "2910"
## [1] "s: " "2911"
## [1] "s: " "2912"
## [1] "s: " "2913"
## [1] "s: " "2914"
## [1] "s: " "2915"
## [1] "s: " "2916"
## [1] "s: " "2917"
## [1] "s: " "2918"
## [1] "s: " "2919"
## [1] "s: " "2920"
## [1] "s: " "2921"
## [1] "s: " "2922"
## [1] "s: " "2923"
## [1] "s: " "2924"
## [1] "s: " "2925"
## [1] "s: " "2926"
## [1] "s: " "2927"
## [1] "s: " "2928"
## [1] "s: " "2929"
## [1] "s: " "2930"
## [1] "s: " "2931"
## [1] "s: " "2932"
## [1] "s: " "2933"
## [1] "s: " "2934"
## [1] "s: " "2935"
## [1] "s: " "2936"
## [1] "s: " "2937"
## [1] "s: " "2938"
## [1] "s: " "2939"
## [1] "s: " "2940"
## [1] "s: " "2941"
## [1] "s: " "2942"
## [1] "s: " "2943"
## [1] "s: " "2944"
## [1] "s: " "2945"
## [1] "s: " "2946"
## [1] "s: " "2947"
## [1] "s: " "2948"
## [1] "s: " "2949"
## [1] "s: " "2950"
## [1] "s: " "2951"
## [1] "s: " "2952"
## [1] "s: " "2953"
## [1] "s: " "2954"
## [1] "s: " "2955"
## [1] "s: " "2956"
## [1] "s: " "2957"
## [1] "s: " "2958"
## [1] "s: " "2959"
## [1] "s: " "2960"
## [1] "s: " "2961"
## [1] "s: " "2962"
## [1] "s: " "2963"
## [1] "s: " "2964"
## [1] "s: " "2965"
## [1] "s: " "2966"
## [1] "s: " "2967"
## [1] "s: " "2968"
## [1] "s: " "2969"
## [1] "s: " "2970"
## [1] "s: " "2971"
## [1] "s: " "2972"
## [1] "s: " "2973"
## [1] "s: " "2974"
## [1] "s: " "2975"
## [1] "s: " "2976"
## [1] "s: " "2977"
## [1] "s: " "2978"
## [1] "s: " "2979"
## [1] "s: " "2980"
## [1] "s: " "2981"
## [1] "s: " "2982"
## [1] "s: " "2983"
## [1] "s: " "2984"
## [1] "s: " "2985"
## [1] "s: " "2986"
## [1] "s: " "2987"
## [1] "s: " "2988"
## [1] "s: " "2989"
## [1] "s: " "2990"
## [1] "s: " "2991"
## [1] "s: " "2992"
## [1] "s: " "2993"
## [1] "s: " "2994"
## [1] "s: " "2995"
## [1] "s: " "2996"
## [1] "s: " "2997"
## [1] "s: " "2998"
## [1] "s: " "2999"
## [1] "s: " "3000"
## [1] "s: " "3001"
## [1] "s: " "3002"
## [1] "s: " "3003"
## [1] "s: " "3004"
## [1] "s: " "3005"
## [1] "s: " "3006"
## [1] "s: " "3007"
## [1] "s: " "3008"
## [1] "s: " "3009"
## [1] "s: " "3010"
## [1] "s: " "3011"
## [1] "s: " "3012"
## [1] "s: " "3013"
## [1] "s: " "3014"
## [1] "s: " "3015"
## [1] "s: " "3016"
## [1] "s: " "3017"
## [1] "s: " "3018"
## [1] "s: " "3019"
## [1] "s: " "3020"
## [1] "s: " "3021"
## [1] "s: " "3022"
## [1] "s: " "3023"
## [1] "s: " "3024"
## [1] "s: " "3025"
## [1] "s: " "3026"
## [1] "s: " "3027"
## [1] "s: " "3028"
## [1] "s: " "3029"
## [1] "s: " "3030"
## [1] "s: " "3031"
## [1] "s: " "3032"
## [1] "s: " "3033"
## [1] "s: " "3034"
## [1] "s: " "3035"
## [1] "s: " "3036"
## [1] "s: " "3037"
## [1] "s: " "3038"
## [1] "s: " "3039"
## [1] "s: " "3040"
## [1] "s: " "3041"
## [1] "s: " "3042"
## [1] "s: " "3043"
## [1] "s: " "3044"
## [1] "s: " "3045"
## [1] "s: " "3046"
## [1] "s: " "3047"
## [1] "s: " "3048"
## [1] "s: " "3049"
## [1] "s: " "3050"
## [1] "s: " "3051"
## [1] "s: " "3052"
## [1] "s: " "3053"
## [1] "s: " "3054"
## [1] "s: " "3055"
## [1] "s: " "3056"
## [1] "s: " "3057"
## [1] "s: " "3058"
## [1] "s: " "3059"
## [1] "s: " "3060"
## [1] "s: " "3061"
## [1] "s: " "3062"
## [1] "s: " "3063"
## [1] "s: " "3064"
## [1] "s: " "3065"
## [1] "s: " "3066"
## [1] "s: " "3067"
## [1] "s: " "3068"
## [1] "s: " "3069"
## [1] "s: " "3070"
## [1] "s: " "3071"
## [1] "s: " "3072"
## [1] "s: " "3073"
## [1] "s: " "3074"
## [1] "s: " "3075"
## [1] "s: " "3076"
## [1] "s: " "3077"
## [1] "s: " "3078"
## [1] "s: " "3079"
## [1] "s: " "3080"
## [1] "s: " "3081"
## [1] "s: " "3082"
## [1] "s: " "3083"
## [1] "s: " "3084"
## [1] "s: " "3085"
## [1] "s: " "3086"
## [1] "s: " "3087"
## [1] "s: " "3088"
## [1] "s: " "3089"
## [1] "s: " "3090"
## [1] "s: " "3091"
## [1] "s: " "3092"
## [1] "s: " "3093"
## [1] "s: " "3094"
## [1] "s: " "3095"
## [1] "s: " "3096"
## [1] "s: " "3097"
## [1] "s: " "3098"
## [1] "s: " "3099"
## [1] "s: " "3100"
## [1] "s: " "3101"
## [1] "s: " "3102"
## [1] "s: " "3103"
## [1] "s: " "3104"
## [1] "s: " "3105"
## [1] "s: " "3106"
## [1] "s: " "3107"
## [1] "s: " "3108"
## [1] "s: " "3109"
## [1] "s: " "3110"
## [1] "s: " "3111"
## [1] "s: " "3112"
## [1] "s: " "3113"
## [1] "s: " "3114"
## [1] "s: " "3115"
## [1] "s: " "3116"
## [1] "s: " "3117"
## [1] "s: " "3118"
## [1] "s: " "3119"
## [1] "s: " "3120"
## [1] "s: " "3121"
## [1] "s: " "3122"
## [1] "s: " "3123"
## [1] "s: " "3124"
## [1] "s: " "3125"
## [1] "s: " "3126"
## [1] "s: " "3127"
## [1] "s: " "3128"
## [1] "s: " "3129"
## [1] "s: " "3130"
## [1] "s: " "3131"
## [1] "s: " "3132"
## [1] "s: " "3133"
## [1] "s: " "3134"
## [1] "s: " "3135"
## [1] "s: " "3136"
## [1] "s: " "3137"
## [1] "s: " "3138"
## [1] "s: " "3139"
## [1] "s: " "3140"
## [1] "s: " "3141"
## [1] "s: " "3142"
## [1] "s: " "3143"
## [1] "s: " "3144"
## [1] "s: " "3145"
## [1] "s: " "3146"
## [1] "s: " "3147"
## [1] "s: " "3148"
## [1] "s: " "3149"
## [1] "s: " "3150"
## [1] "s: " "3151"
## [1] "s: " "3152"
## [1] "s: " "3153"
## [1] "s: " "3154"
## [1] "s: " "3155"
## [1] "s: " "3156"
## [1] "s: " "3157"
## [1] "s: " "3158"
## [1] "s: " "3159"
## [1] "s: " "3160"
## [1] "s: " "3161"
## [1] "s: " "3162"
## [1] "s: " "3163"
## [1] "s: " "3164"
## [1] "s: " "3165"
## [1] "s: " "3166"
## [1] "s: " "3167"
## [1] "s: " "3168"
## [1] "s: " "3169"
## [1] "s: " "3170"
## [1] "s: " "3171"
## [1] "s: " "3172"
## [1] "s: " "3173"
## [1] "s: " "3174"
## [1] "s: " "3175"
## [1] "s: " "3176"
## [1] "s: " "3177"
## [1] "s: " "3178"
## [1] "s: " "3179"
## [1] "s: " "3180"
## [1] "s: " "3181"
## [1] "s: " "3182"
## [1] "s: " "3183"
## [1] "s: " "3184"
## [1] "s: " "3185"
## [1] "s: " "3186"
## [1] "s: " "3187"
## [1] "s: " "3188"
## [1] "s: " "3189"
## [1] "s: " "3190"
## [1] "s: " "3191"
## [1] "s: " "3192"
## [1] "s: " "3193"
## [1] "s: " "3194"
## [1] "s: " "3195"
## [1] "s: " "3196"
## [1] "s: " "3197"
## [1] "s: " "3198"
## [1] "s: " "3199"
## [1] "s: " "3200"
## [1] "s: " "3201"
## [1] "s: " "3202"
## [1] "s: " "3203"
## [1] "s: " "3204"
## [1] "s: " "3205"
## [1] "s: " "3206"
## [1] "s: " "3207"
## [1] "s: " "3208"
## [1] "s: " "3209"
## [1] "s: " "3210"
## [1] "s: " "3211"
## [1] "s: " "3212"
## [1] "s: " "3213"
## [1] "s: " "3214"
## [1] "s: " "3215"
## [1] "s: " "3216"
## [1] "s: " "3217"
## [1] "s: " "3218"
## [1] "s: " "3219"
## [1] "s: " "3220"
## [1] "s: " "3221"
## [1] "s: " "3222"
## [1] "s: " "3223"
## [1] "s: " "3224"
## [1] "s: " "3225"
## [1] "s: " "3226"
## [1] "s: " "3227"
## [1] "s: " "3228"
## [1] "s: " "3229"
## [1] "s: " "3230"
## [1] "s: " "3231"
## [1] "s: " "3232"
## [1] "s: " "3233"
## [1] "s: " "3234"
## [1] "s: " "3235"
## [1] "s: " "3236"
## [1] "s: " "3237"
## [1] "s: " "3238"
## [1] "s: " "3239"
## [1] "s: " "3240"
## [1] "s: " "3241"
## [1] "s: " "3242"
## [1] "s: " "3243"
## [1] "s: " "3244"
## [1] "s: " "3245"
## [1] "s: " "3246"
## [1] "s: " "3247"
## [1] "s: " "3248"
## [1] "s: " "3249"
## [1] "s: " "3250"
## [1] "s: " "3251"
## [1] "s: " "3252"
## [1] "s: " "3253"
## [1] "s: " "3254"
## [1] "s: " "3255"
## [1] "s: " "3256"
## [1] "s: " "3257"
## [1] "s: " "3258"
## [1] "s: " "3259"
## [1] "s: " "3260"
## [1] "s: " "3261"
## [1] "s: " "3262"
## [1] "s: " "3263"
## [1] "s: " "3264"
## [1] "s: " "3265"
## [1] "s: " "3266"
## [1] "s: " "3267"
## [1] "s: " "3268"
## [1] "s: " "3269"
## [1] "s: " "3270"
## [1] "s: " "3271"
## [1] "s: " "3272"
## [1] "s: " "3273"
## [1] "s: " "3274"
## [1] "s: " "3275"
## [1] "s: " "3276"
## [1] "s: " "3277"
## [1] "s: " "3278"
## [1] "s: " "3279"
## [1] "s: " "3280"
## [1] "s: " "3281"
## [1] "s: " "3282"
## [1] "s: " "3283"
## [1] "s: " "3284"
## [1] "s: " "3285"
## [1] "s: " "3286"
## [1] "s: " "3287"
## [1] "s: " "3288"
## [1] "s: " "3289"
## [1] "s: " "3290"
## [1] "s: " "3291"
## [1] "s: " "3292"
## [1] "s: " "3293"
## [1] "s: " "3294"
## [1] "s: " "3295"
## [1] "s: " "3296"
## [1] "s: " "3297"
## [1] "s: " "3298"
## [1] "s: " "3299"
## [1] "s: " "3300"
## [1] "s: " "3301"
## [1] "s: " "3302"
## [1] "s: " "3303"
## [1] "s: " "3304"
## [1] "s: " "3305"
## [1] "s: " "3306"
## [1] "s: " "3307"
## [1] "s: " "3308"
## [1] "s: " "3309"
## [1] "s: " "3310"
## [1] "s: " "3311"
## [1] "s: " "3312"
## [1] "s: " "3313"
## [1] "s: " "3314"
## [1] "s: " "3315"
## [1] "s: " "3316"
## [1] "s: " "3317"
## [1] "s: " "3318"
## [1] "s: " "3319"
## [1] "s: " "3320"
## [1] "s: " "3321"
## [1] "s: " "3322"
## [1] "s: " "3323"
## [1] "s: " "3324"
## [1] "s: " "3325"
## [1] "s: " "3326"
## [1] "s: " "3327"
## [1] "s: " "3328"
## [1] "s: " "3329"
## [1] "s: " "3330"
## [1] "s: " "3331"
## [1] "s: " "3332"
## [1] "s: " "3333"
## [1] "s: " "3334"
## [1] "s: " "3335"
## [1] "s: " "3336"
## [1] "s: " "3337"
## [1] "s: " "3338"
## [1] "s: " "3339"
## [1] "s: " "3340"
## [1] "s: " "3341"
## [1] "s: " "3342"
## [1] "s: " "3343"
## [1] "s: " "3344"
## [1] "s: " "3345"
## [1] "s: " "3346"
## [1] "s: " "3347"
## [1] "s: " "3348"
## [1] "s: " "3349"
## [1] "s: " "3350"
## [1] "s: " "3351"
## [1] "s: " "3352"
## [1] "s: " "3353"
## [1] "s: " "3354"
## [1] "s: " "3355"
## [1] "s: " "3356"
## [1] "s: " "3357"
## [1] "s: " "3358"
## [1] "s: " "3359"
## [1] "s: " "3360"
## [1] "s: " "3361"
## [1] "s: " "3362"
## [1] "s: " "3363"
## [1] "s: " "3364"
## [1] "s: " "3365"
## [1] "s: " "3366"
## [1] "s: " "3367"
## [1] "s: " "3368"
## [1] "s: " "3369"
## [1] "s: " "3370"
## [1] "s: " "3371"
## [1] "s: " "3372"
## [1] "s: " "3373"
## [1] "s: " "3374"
## [1] "s: " "3375"
## [1] "s: " "3376"
## [1] "s: " "3377"
## [1] "s: " "3378"
## [1] "s: " "3379"
## [1] "s: " "3380"
## [1] "s: " "3381"
## [1] "s: " "3382"
## [1] "s: " "3383"
## [1] "s: " "3384"
## [1] "s: " "3385"
## [1] "s: " "3386"
## [1] "s: " "3387"
## [1] "s: " "3388"
## [1] "s: " "3389"
## [1] "s: " "3390"
## [1] "s: " "3391"
## [1] "s: " "3392"
## [1] "s: " "3393"
## [1] "s: " "3394"
## [1] "s: " "3395"
## [1] "s: " "3396"
## [1] "s: " "3397"
## [1] "s: " "3398"
## [1] "s: " "3399"
## [1] "s: " "3400"
## [1] "s: " "3401"
## [1] "s: " "3402"
## [1] "s: " "3403"
## [1] "s: " "3404"
## [1] "s: " "3405"
## [1] "s: " "3406"
## [1] "s: " "3407"
## [1] "s: " "3408"
## [1] "s: " "3409"
## [1] "s: " "3410"
## [1] "s: " "3411"
## [1] "s: " "3412"
## [1] "s: " "3413"
## [1] "s: " "3414"
## [1] "s: " "3415"
## [1] "s: " "3416"
## [1] "s: " "3417"
## [1] "s: " "3418"
## [1] "s: " "3419"
## [1] "s: " "3420"
## [1] "s: " "3421"
## [1] "s: " "3422"
## [1] "s: " "3423"
## [1] "s: " "3424"
## [1] "s: " "3425"
## [1] "s: " "3426"
## [1] "s: " "3427"
## [1] "s: " "3428"
## [1] "s: " "3429"
## [1] "s: " "3430"
## [1] "s: " "3431"
## [1] "s: " "3432"
## [1] "s: " "3433"
## [1] "s: " "3434"
## [1] "s: " "3435"
## [1] "s: " "3436"
## [1] "s: " "3437"
## [1] "s: " "3438"
## [1] "s: " "3439"
## [1] "s: " "3440"
## [1] "s: " "3441"
## [1] "s: " "3442"
## [1] "s: " "3443"
## [1] "s: " "3444"
## [1] "s: " "3445"
## [1] "s: " "3446"
## [1] "s: " "3447"
## [1] "s: " "3448"
## [1] "s: " "3449"
## [1] "s: " "3450"
## [1] "s: " "3451"
## [1] "s: " "3452"
## [1] "s: " "3453"
## [1] "s: " "3454"
## [1] "s: " "3455"
## [1] "s: " "3456"
## [1] "s: " "3457"
## [1] "s: " "3458"
## [1] "s: " "3459"
## [1] "s: " "3460"
## [1] "s: " "3461"
## [1] "s: " "3462"
## [1] "s: " "3463"
## [1] "s: " "3464"
## [1] "s: " "3465"
## [1] "s: " "3466"
## [1] "s: " "3467"
## [1] "s: " "3468"
## [1] "s: " "3469"
## [1] "s: " "3470"
## [1] "s: " "3471"
## [1] "s: " "3472"
## [1] "s: " "3473"
## [1] "s: " "3474"
## [1] "s: " "3475"
## [1] "s: " "3476"
## [1] "s: " "3477"
## [1] "s: " "3478"
## [1] "s: " "3479"
## [1] "s: " "3480"
## [1] "s: " "3481"
## [1] "s: " "3482"
## [1] "s: " "3483"
## [1] "s: " "3484"
## [1] "s: " "3485"
## [1] "s: " "3486"
## [1] "s: " "3487"
## [1] "s: " "3488"
## [1] "s: " "3489"
## [1] "s: " "3490"
## [1] "s: " "3491"
## [1] "s: " "3492"
## [1] "s: " "3493"
## [1] "s: " "3494"
## [1] "s: " "3495"
## [1] "s: " "3496"
## [1] "s: " "3497"
## [1] "s: " "3498"
## [1] "s: " "3499"
## [1] "s: " "3500"
## [1] "s: " "3501"
## [1] "s: " "3502"
## [1] "s: " "3503"
## [1] "s: " "3504"
## [1] "s: " "3505"
## [1] "s: " "3506"
## [1] "s: " "3507"
## [1] "s: " "3508"
## [1] "s: " "3509"
## [1] "s: " "3510"
## [1] "s: " "3511"
## [1] "s: " "3512"
## [1] "s: " "3513"
## [1] "s: " "3514"
## [1] "s: " "3515"
## [1] "s: " "3516"
## [1] "s: " "3517"
## [1] "s: " "3518"
## [1] "s: " "3519"
## [1] "s: " "3520"
## [1] "s: " "3521"
## [1] "s: " "3522"
## [1] "s: " "3523"
## [1] "s: " "3524"
## [1] "s: " "3525"
## [1] "s: " "3526"
## [1] "s: " "3527"
## [1] "s: " "3528"
## [1] "s: " "3529"
## [1] "s: " "3530"
## [1] "s: " "3531"
## [1] "s: " "3532"
## [1] "s: " "3533"
## [1] "s: " "3534"
## [1] "s: " "3535"
## [1] "s: " "3536"
## [1] "s: " "3537"
## [1] "s: " "3538"
## [1] "s: " "3539"
## [1] "s: " "3540"
## [1] "s: " "3541"
## [1] "s: " "3542"
## [1] "s: " "3543"
## [1] "s: " "3544"
## [1] "s: " "3545"
## [1] "s: " "3546"
## [1] "s: " "3547"
## [1] "s: " "3548"
## [1] "s: " "3549"
## [1] "s: " "3550"
## [1] "s: " "3551"
## [1] "s: " "3552"
## [1] "s: " "3553"
## [1] "s: " "3554"
## [1] "s: " "3555"
## [1] "s: " "3556"
## [1] "s: " "3557"
## [1] "s: " "3558"
## [1] "s: " "3559"
## [1] "s: " "3560"
## [1] "s: " "3561"
## [1] "s: " "3562"
## [1] "s: " "3563"
## [1] "s: " "3564"
## [1] "s: " "3565"
## [1] "s: " "3566"
## [1] "s: " "3567"
## [1] "s: " "3568"
## [1] "s: " "3569"
## [1] "s: " "3570"
## [1] "s: " "3571"
## [1] "s: " "3572"
## [1] "s: " "3573"
## [1] "s: " "3574"
## [1] "s: " "3575"
## [1] "s: " "3576"
## [1] "s: " "3577"
## [1] "s: " "3578"
## [1] "s: " "3579"
## [1] "s: " "3580"
## [1] "s: " "3581"
## [1] "s: " "3582"
## [1] "s: " "3583"
## [1] "s: " "3584"
## [1] "s: " "3585"
## [1] "s: " "3586"
## [1] "s: " "3587"
## [1] "s: " "3588"
## [1] "s: " "3589"
## [1] "s: " "3590"
## [1] "s: " "3591"
## [1] "s: " "3592"
## [1] "s: " "3593"
## [1] "s: " "3594"
## [1] "s: " "3595"
## [1] "s: " "3596"
## [1] "s: " "3597"
## [1] "s: " "3598"
## [1] "s: " "3599"
## [1] "s: " "3600"
## [1] "s: " "3601"
## [1] "s: " "3602"
## [1] "s: " "3603"
## [1] "s: " "3604"
## [1] "s: " "3605"
## [1] "s: " "3606"
## [1] "s: " "3607"
## [1] "s: " "3608"
## [1] "s: " "3609"
## [1] "s: " "3610"
## [1] "s: " "3611"
## [1] "s: " "3612"
## [1] "s: " "3613"
## [1] "s: " "3614"
## [1] "s: " "3615"
## [1] "s: " "3616"
## [1] "s: " "3617"
## [1] "s: " "3618"
## [1] "s: " "3619"
## [1] "s: " "3620"
## [1] "s: " "3621"
## [1] "s: " "3622"
## [1] "s: " "3623"
## [1] "s: " "3624"
## [1] "s: " "3625"
## [1] "s: " "3626"
## [1] "s: " "3627"
## [1] "s: " "3628"
## [1] "s: " "3629"
## [1] "s: " "3630"
## [1] "s: " "3631"
## [1] "s: " "3632"
## [1] "s: " "3633"
## [1] "s: " "3634"
## [1] "s: " "3635"
## [1] "s: " "3636"
## [1] "s: " "3637"
## [1] "s: " "3638"
## [1] "s: " "3639"
## [1] "s: " "3640"
## [1] "s: " "3641"
## [1] "s: " "3642"
## [1] "s: " "3643"
## [1] "s: " "3644"
## [1] "s: " "3645"
## [1] "s: " "3646"
## [1] "s: " "3647"
## [1] "s: " "3648"
## [1] "s: " "3649"
## [1] "s: " "3650"
## [1] "s: " "3651"
## [1] "s: " "3652"
## [1] "s: " "3653"
## [1] "s: " "3654"
## [1] "s: " "3655"
## [1] "s: " "3656"
## [1] "s: " "3657"
## [1] "s: " "3658"
## [1] "s: " "3659"
## [1] "s: " "3660"
## [1] "s: " "3661"
## [1] "s: " "3662"
## [1] "s: " "3663"
## [1] "s: " "3664"
## [1] "s: " "3665"
## [1] "s: " "3666"
## [1] "s: " "3667"
## [1] "s: " "3668"
## [1] "s: " "3669"
## [1] "s: " "3670"
## [1] "s: " "3671"
## [1] "s: " "3672"
## [1] "s: " "3673"
## [1] "s: " "3674"
## [1] "s: " "3675"
## [1] "s: " "3676"
## [1] "s: " "3677"
## [1] "s: " "3678"
## [1] "s: " "3679"
## [1] "s: " "3680"
## [1] "s: " "3681"
## [1] "s: " "3682"
## [1] "s: " "3683"
## [1] "s: " "3684"
## [1] "s: " "3685"
## [1] "s: " "3686"
## [1] "s: " "3687"
## [1] "s: " "3688"
## [1] "s: " "3689"
## [1] "s: " "3690"
## [1] "s: " "3691"
## [1] "s: " "3692"
## [1] "s: " "3693"
## [1] "s: " "3694"
## [1] "s: " "3695"
## [1] "s: " "3696"
## [1] "s: " "3697"
## [1] "s: " "3698"
## [1] "s: " "3699"
## [1] "s: " "3700"
## [1] "s: " "3701"
## [1] "s: " "3702"
## [1] "s: " "3703"
## [1] "s: " "3704"
## [1] "s: " "3705"
## [1] "s: " "3706"
## [1] "s: " "3707"
## [1] "s: " "3708"
## [1] "s: " "3709"
## [1] "s: " "3710"
## [1] "s: " "3711"
## [1] "s: " "3712"
## [1] "s: " "3713"
## [1] "s: " "3714"
## [1] "s: " "3715"
## [1] "s: " "3716"
## [1] "s: " "3717"
## [1] "s: " "3718"
## [1] "s: " "3719"
## [1] "s: " "3720"
## [1] "s: " "3721"
## [1] "s: " "3722"
## [1] "s: " "3723"
## [1] "s: " "3724"
## [1] "s: " "3725"
## [1] "s: " "3726"
## [1] "s: " "3727"
## [1] "s: " "3728"
## [1] "s: " "3729"
## [1] "s: " "3730"
## [1] "s: " "3731"
## [1] "s: " "3732"
## [1] "s: " "3733"
## [1] "s: " "3734"
## [1] "s: " "3735"
## [1] "s: " "3736"
## [1] "s: " "3737"
## [1] "s: " "3738"
## [1] "s: " "3739"
## [1] "s: " "3740"
## [1] "s: " "3741"
## [1] "s: " "3742"
## [1] "s: " "3743"
## [1] "s: " "3744"
## [1] "s: " "3745"
## [1] "s: " "3746"
## [1] "s: " "3747"
## [1] "s: " "3748"
## [1] "s: " "3749"
## [1] "s: " "3750"
## [1] "s: " "3751"
## [1] "s: " "3752"
## [1] "s: " "3753"
## [1] "s: " "3754"
## [1] "s: " "3755"
## [1] "s: " "3756"
## [1] "s: " "3757"
## [1] "s: " "3758"
## [1] "s: " "3759"
## [1] "s: " "3760"
## [1] "s: " "3761"
## [1] "s: " "3762"
## [1] "s: " "3763"
## [1] "s: " "3764"
## [1] "s: " "3765"
## [1] "s: " "3766"
## [1] "s: " "3767"
## [1] "s: " "3768"
## [1] "s: " "3769"
## [1] "s: " "3770"
## [1] "s: " "3771"
## [1] "s: " "3772"
## [1] "s: " "3773"
## [1] "s: " "3774"
## [1] "s: " "3775"
## [1] "s: " "3776"
## [1] "s: " "3777"
## [1] "s: " "3778"
## [1] "s: " "3779"
## [1] "s: " "3780"
## [1] "s: " "3781"
## [1] "s: " "3782"
## [1] "s: " "3783"
## [1] "s: " "3784"
## [1] "s: " "3785"
## [1] "s: " "3786"
## [1] "s: " "3787"
## [1] "s: " "3788"
## [1] "s: " "3789"
## [1] "s: " "3790"
## [1] "s: " "3791"
## [1] "s: " "3792"
## [1] "s: " "3793"
## [1] "s: " "3794"
## [1] "s: " "3795"
## [1] "s: " "3796"
## [1] "s: " "3797"
## [1] "s: " "3798"
## [1] "s: " "3799"
## [1] "s: " "3800"
## [1] "s: " "3801"
## [1] "s: " "3802"
## [1] "s: " "3803"
## [1] "s: " "3804"
## [1] "s: " "3805"
## [1] "s: " "3806"
## [1] "s: " "3807"
## [1] "s: " "3808"
## [1] "s: " "3809"
## [1] "s: " "3810"
## [1] "s: " "3811"
## [1] "s: " "3812"
## [1] "s: " "3813"
## [1] "s: " "3814"
## [1] "s: " "3815"
## [1] "s: " "3816"
## [1] "s: " "3817"
## [1] "s: " "3818"
## [1] "s: " "3819"
## [1] "s: " "3820"
## [1] "s: " "3821"
## [1] "s: " "3822"
## [1] "s: " "3823"
## [1] "s: " "3824"
## [1] "s: " "3825"
## [1] "s: " "3826"
## [1] "s: " "3827"
## [1] "s: " "3828"
## [1] "s: " "3829"
## [1] "s: " "3830"
## [1] "s: " "3831"
## [1] "s: " "3832"
## [1] "s: " "3833"
## [1] "s: " "3834"
## [1] "s: " "3835"
## [1] "s: " "3836"
## [1] "s: " "3837"
## [1] "s: " "3838"
## [1] "s: " "3839"
## [1] "s: " "3840"
## [1] "s: " "3841"
## [1] "s: " "3842"
## [1] "s: " "3843"
## [1] "s: " "3844"
## [1] "s: " "3845"
## [1] "s: " "3846"
## [1] "s: " "3847"
## [1] "s: " "3848"
## [1] "s: " "3849"
## [1] "s: " "3850"
## [1] "s: " "3851"
## [1] "s: " "3852"
## [1] "s: " "3853"
## [1] "s: " "3854"
## [1] "s: " "3855"
## [1] "s: " "3856"
## [1] "s: " "3857"
## [1] "s: " "3858"
## [1] "s: " "3859"
## [1] "s: " "3860"
## [1] "s: " "3861"
## [1] "s: " "3862"
## [1] "s: " "3863"
## [1] "s: " "3864"
## [1] "s: " "3865"
## [1] "s: " "3866"
## [1] "s: " "3867"
## [1] "s: " "3868"
## [1] "s: " "3869"
## [1] "s: " "3870"
## [1] "s: " "3871"
## [1] "s: " "3872"
## [1] "s: " "3873"
## [1] "s: " "3874"
## [1] "s: " "3875"
## [1] "s: " "3876"
## [1] "s: " "3877"
## [1] "s: " "3878"
## [1] "s: " "3879"
## [1] "s: " "3880"
## [1] "s: " "3881"
## [1] "s: " "3882"
## [1] "s: " "3883"
## [1] "s: " "3884"
## [1] "s: " "3885"
## [1] "s: " "3886"
## [1] "s: " "3887"
## [1] "s: " "3888"
## [1] "s: " "3889"
## [1] "s: " "3890"
## [1] "s: " "3891"
## [1] "s: " "3892"
## [1] "s: " "3893"
## [1] "s: " "3894"
## [1] "s: " "3895"
## [1] "s: " "3896"
## [1] "s: " "3897"
## [1] "s: " "3898"
## [1] "s: " "3899"
## [1] "s: " "3900"
## [1] "s: " "3901"
## [1] "s: " "3902"
## [1] "s: " "3903"
## [1] "s: " "3904"
## [1] "s: " "3905"
## [1] "s: " "3906"
## [1] "s: " "3907"
## [1] "s: " "3908"
## [1] "s: " "3909"
## [1] "s: " "3910"
## [1] "s: " "3911"
## [1] "s: " "3912"
## [1] "s: " "3913"
## [1] "s: " "3914"
## [1] "s: " "3915"
## [1] "s: " "3916"
## [1] "s: " "3917"
## [1] "s: " "3918"
## [1] "s: " "3919"
## [1] "s: " "3920"
## [1] "s: " "3921"
## [1] "s: " "3922"
## [1] "s: " "3923"
## [1] "s: " "3924"
## [1] "s: " "3925"
## [1] "s: " "3926"
## [1] "s: " "3927"
## [1] "s: " "3928"
## [1] "s: " "3929"
## [1] "s: " "3930"
## [1] "s: " "3931"
## [1] "s: " "3932"
## [1] "s: " "3933"
## [1] "s: " "3934"
## [1] "s: " "3935"
## [1] "s: " "3936"
## [1] "s: " "3937"
## [1] "s: " "3938"
## [1] "s: " "3939"
## [1] "s: " "3940"
## [1] "s: " "3941"
## [1] "s: " "3942"
## [1] "s: " "3943"
## [1] "s: " "3944"
## [1] "s: " "3945"
## [1] "s: " "3946"
## [1] "s: " "3947"
## [1] "s: " "3948"
## [1] "s: " "3949"
## [1] "s: " "3950"
## [1] "s: " "3951"
## [1] "s: " "3952"
## [1] "s: " "3953"
## [1] "s: " "3954"
## [1] "s: " "3955"
## [1] "s: " "3956"
## [1] "s: " "3957"
## [1] "s: " "3958"
## [1] "s: " "3959"
## [1] "s: " "3960"
## [1] "s: " "3961"
## [1] "s: " "3962"
## [1] "s: " "3963"
## [1] "s: " "3964"
## [1] "s: " "3965"
## [1] "s: " "3966"
## [1] "s: " "3967"
## [1] "s: " "3968"
## [1] "s: " "3969"
## [1] "s: " "3970"
## [1] "s: " "3971"
## [1] "s: " "3972"
## [1] "s: " "3973"
## [1] "s: " "3974"
## [1] "s: " "3975"
## [1] "s: " "3976"
## [1] "s: " "3977"
## [1] "s: " "3978"
## [1] "s: " "3979"
## [1] "s: " "3980"
## [1] "s: " "3981"
## [1] "s: " "3982"
## [1] "s: " "3983"
## [1] "s: " "3984"
## [1] "s: " "3985"
## [1] "s: " "3986"
## [1] "s: " "3987"
## [1] "s: " "3988"
## [1] "s: " "3989"
## [1] "s: " "3990"
## [1] "s: " "3991"
## [1] "s: " "3992"
## [1] "s: " "3993"
## [1] "s: " "3994"
## [1] "s: " "3995"
## [1] "s: " "3996"
## [1] "s: " "3997"
## [1] "s: " "3998"
## [1] "s: " "3999"
## [1] "s: " "4000"
## [1] "s: " "4001"
## [1] "s: " "4002"
## [1] "s: " "4003"
## [1] "s: " "4004"
## [1] "s: " "4005"
## [1] "s: " "4006"
## [1] "s: " "4007"
## [1] "s: " "4008"
## [1] "s: " "4009"
## [1] "s: " "4010"
## [1] "s: " "4011"
## [1] "s: " "4012"
## [1] "s: " "4013"
## [1] "s: " "4014"
## [1] "s: " "4015"
## [1] "s: " "4016"
## [1] "s: " "4017"
## [1] "s: " "4018"
## [1] "s: " "4019"
## [1] "s: " "4020"
## [1] "s: " "4021"
## [1] "s: " "4022"
## [1] "s: " "4023"
## [1] "s: " "4024"
## [1] "s: " "4025"
## [1] "s: " "4026"
## [1] "s: " "4027"
## [1] "s: " "4028"
## [1] "s: " "4029"
## [1] "s: " "4030"
## [1] "s: " "4031"
## [1] "s: " "4032"
## [1] "s: " "4033"
## [1] "s: " "4034"
## [1] "s: " "4035"
## [1] "s: " "4036"
## [1] "s: " "4037"
## [1] "s: " "4038"
## [1] "s: " "4039"
## [1] "s: " "4040"
## [1] "s: " "4041"
## [1] "s: " "4042"
## [1] "s: " "4043"
## [1] "s: " "4044"
## [1] "s: " "4045"
## [1] "s: " "4046"
## [1] "s: " "4047"
## [1] "s: " "4048"
## [1] "s: " "4049"
## [1] "s: " "4050"
## [1] "s: " "4051"
## [1] "s: " "4052"
## [1] "s: " "4053"
## [1] "s: " "4054"
## [1] "s: " "4055"
## [1] "s: " "4056"
## [1] "s: " "4057"
## [1] "s: " "4058"
## [1] "s: " "4059"
## [1] "s: " "4060"
## [1] "s: " "4061"
## [1] "s: " "4062"
## [1] "s: " "4063"
## [1] "s: " "4064"
## [1] "s: " "4065"
## [1] "s: " "4066"
## [1] "s: " "4067"
## [1] "s: " "4068"
## [1] "s: " "4069"
## [1] "s: " "4070"
## [1] "s: " "4071"
## [1] "s: " "4072"
## [1] "s: " "4073"
## [1] "s: " "4074"
## [1] "s: " "4075"
## [1] "s: " "4076"
## [1] "s: " "4077"
## [1] "s: " "4078"
## [1] "s: " "4079"
## [1] "s: " "4080"
## [1] "s: " "4081"
## [1] "s: " "4082"
## [1] "s: " "4083"
## [1] "s: " "4084"
## [1] "s: " "4085"
## [1] "s: " "4086"
## [1] "s: " "4087"
## [1] "s: " "4088"
## [1] "s: " "4089"
## [1] "s: " "4090"
## [1] "s: " "4091"
## [1] "s: " "4092"
## [1] "s: " "4093"
## [1] "s: " "4094"
## [1] "s: " "4095"
## [1] "s: " "4096"
## [1] "s: " "4097"
## [1] "s: " "4098"
## [1] "s: " "4099"
## [1] "s: " "4100"
## [1] "s: " "4101"
## [1] "s: " "4102"
## [1] "s: " "4103"
## [1] "s: " "4104"
## [1] "s: " "4105"
## [1] "s: " "4106"
## [1] "s: " "4107"
## [1] "s: " "4108"
## [1] "s: " "4109"
## [1] "s: " "4110"
## [1] "s: " "4111"
## [1] "s: " "4112"
## [1] "s: " "4113"
## [1] "s: " "4114"
## [1] "s: " "4115"
## [1] "s: " "4116"
## [1] "s: " "4117"
## [1] "s: " "4118"
## [1] "s: " "4119"
## [1] "s: " "4120"
## [1] "s: " "4121"
## [1] "s: " "4122"
## [1] "s: " "4123"
## [1] "s: " "4124"
## [1] "s: " "4125"
## [1] "s: " "4126"
## [1] "s: " "4127"
## [1] "s: " "4128"
## [1] "s: " "4129"
## [1] "s: " "4130"
## [1] "s: " "4131"
## [1] "s: " "4132"
## [1] "s: " "4133"
## [1] "s: " "4134"
## [1] "s: " "4135"
## [1] "s: " "4136"
## [1] "s: " "4137"
## [1] "s: " "4138"
## [1] "s: " "4139"
## [1] "s: " "4140"
## [1] "s: " "4141"
## [1] "s: " "4142"
## [1] "s: " "4143"
## [1] "s: " "4144"
## [1] "s: " "4145"
## [1] "s: " "4146"
## [1] "s: " "4147"
## [1] "s: " "4148"
## [1] "s: " "4149"
## [1] "s: " "4150"
## [1] "s: " "4151"
## [1] "s: " "4152"
## [1] "s: " "4153"
## [1] "s: " "4154"
## [1] "s: " "4155"
## [1] "s: " "4156"
## [1] "s: " "4157"
## [1] "s: " "4158"
## [1] "s: " "4159"
## [1] "s: " "4160"
## [1] "s: " "4161"
## [1] "s: " "4162"
## [1] "s: " "4163"
## [1] "s: " "4164"
## [1] "s: " "4165"
## [1] "s: " "4166"
## [1] "s: " "4167"
## [1] "s: " "4168"
## [1] "s: " "4169"
## [1] "s: " "4170"
## [1] "s: " "4171"
## [1] "s: " "4172"
## [1] "s: " "4173"
## [1] "s: " "4174"
## [1] "s: " "4175"
## [1] "s: " "4176"
## [1] "s: " "4177"
## [1] "s: " "4178"
## [1] "s: " "4179"
## [1] "s: " "4180"
## [1] "s: " "4181"
## [1] "s: " "4182"
## [1] "s: " "4183"
## [1] "s: " "4184"
## [1] "s: " "4185"
## [1] "s: " "4186"
## [1] "s: " "4187"
## [1] "s: " "4188"
## [1] "s: " "4189"
## [1] "s: " "4190"
## [1] "s: " "4191"
## [1] "s: " "4192"
## [1] "s: " "4193"
## [1] "s: " "4194"
## [1] "s: " "4195"
## [1] "s: " "4196"
## [1] "s: " "4197"
## [1] "s: " "4198"
## [1] "s: " "4199"
## [1] "s: " "4200"
## [1] "s: " "4201"
## [1] "s: " "4202"
## [1] "s: " "4203"
## [1] "s: " "4204"
## [1] "s: " "4205"
## [1] "s: " "4206"
## [1] "s: " "4207"
## [1] "s: " "4208"
## [1] "s: " "4209"
## [1] "s: " "4210"
## [1] "s: " "4211"
## [1] "s: " "4212"
## [1] "s: " "4213"
## [1] "s: " "4214"
## [1] "s: " "4215"
## [1] "s: " "4216"
## [1] "s: " "4217"
## [1] "s: " "4218"
## [1] "s: " "4219"
## [1] "s: " "4220"
## [1] "s: " "4221"
## [1] "s: " "4222"
## [1] "s: " "4223"
## [1] "s: " "4224"
## [1] "s: " "4225"
## [1] "s: " "4226"
## [1] "s: " "4227"
## [1] "s: " "4228"
## [1] "s: " "4229"
## [1] "s: " "4230"
## [1] "s: " "4231"
## [1] "s: " "4232"
## [1] "s: " "4233"
## [1] "s: " "4234"
## [1] "s: " "4235"
## [1] "s: " "4236"
## [1] "s: " "4237"
## [1] "s: " "4238"
## [1] "s: " "4239"
## [1] "s: " "4240"
## [1] "s: " "4241"
## [1] "s: " "4242"
## [1] "s: " "4243"
## [1] "s: " "4244"
## [1] "s: " "4245"
## [1] "s: " "4246"
## [1] "s: " "4247"
## [1] "s: " "4248"
## [1] "s: " "4249"
## [1] "s: " "4250"
## [1] "s: " "4251"
## [1] "s: " "4252"
## [1] "s: " "4253"
## [1] "s: " "4254"
## [1] "s: " "4255"
## [1] "s: " "4256"
## [1] "s: " "4257"
## [1] "s: " "4258"
## [1] "s: " "4259"
## [1] "s: " "4260"
## [1] "s: " "4261"
## [1] "s: " "4262"
## [1] "s: " "4263"
## [1] "s: " "4264"
## [1] "s: " "4265"
## [1] "s: " "4266"
## [1] "s: " "4267"
## [1] "s: " "4268"
## [1] "s: " "4269"
## [1] "s: " "4270"
## [1] "s: " "4271"
## [1] "s: " "4272"
## [1] "s: " "4273"
## [1] "s: " "4274"
## [1] "s: " "4275"
## [1] "s: " "4276"
## [1] "s: " "4277"
## [1] "s: " "4278"
## [1] "s: " "4279"
## [1] "s: " "4280"
## [1] "s: " "4281"
## [1] "s: " "4282"
## [1] "s: " "4283"
## [1] "s: " "4284"
## [1] "s: " "4285"
## [1] "s: " "4286"
## [1] "s: " "4287"
## [1] "s: " "4288"
## [1] "s: " "4289"
## [1] "s: " "4290"
## [1] "s: " "4291"
## [1] "s: " "4292"
## [1] "s: " "4293"
## [1] "s: " "4294"
## [1] "s: " "4295"
## [1] "s: " "4296"
## [1] "s: " "4297"
## [1] "s: " "4298"
## [1] "s: " "4299"
## [1] "s: " "4300"
## [1] "s: " "4301"
## [1] "s: " "4302"
## [1] "s: " "4303"
## [1] "s: " "4304"
## [1] "s: " "4305"
## [1] "s: " "4306"
## [1] "s: " "4307"
## [1] "s: " "4308"
## [1] "s: " "4309"
## [1] "s: " "4310"
## [1] "s: " "4311"
## [1] "s: " "4312"
## [1] "s: " "4313"
## [1] "s: " "4314"
## [1] "s: " "4315"
## [1] "s: " "4316"
## [1] "s: " "4317"
## [1] "s: " "4318"
## [1] "s: " "4319"
## [1] "s: " "4320"
## [1] "s: " "4321"
## [1] "s: " "4322"
## [1] "s: " "4323"
## [1] "s: " "4324"
## [1] "s: " "4325"
## [1] "s: " "4326"
## [1] "s: " "4327"
## [1] "s: " "4328"
## [1] "s: " "4329"
## [1] "s: " "4330"
## [1] "s: " "4331"
## [1] "s: " "4332"
## [1] "s: " "4333"
## [1] "s: " "4334"
## [1] "s: " "4335"
## [1] "s: " "4336"
## [1] "s: " "4337"
## [1] "s: " "4338"
## [1] "s: " "4339"
## [1] "s: " "4340"
## [1] "s: " "4341"
## [1] "s: " "4342"
## [1] "s: " "4343"
## [1] "s: " "4344"
## [1] "s: " "4345"
## [1] "s: " "4346"
## [1] "s: " "4347"
## [1] "s: " "4348"
## [1] "s: " "4349"
## [1] "s: " "4350"
## [1] "s: " "4351"
## [1] "s: " "4352"
## [1] "s: " "4353"
## [1] "s: " "4354"
## [1] "s: " "4355"
## [1] "s: " "4356"
## [1] "s: " "4357"
## [1] "s: " "4358"
## [1] "s: " "4359"
## [1] "s: " "4360"
## [1] "s: " "4361"
## [1] "s: " "4362"
## [1] "s: " "4363"
## [1] "s: " "4364"
## [1] "s: " "4365"
## [1] "s: " "4366"
## [1] "s: " "4367"
## [1] "s: " "4368"
## [1] "s: " "4369"
## [1] "s: " "4370"
## [1] "s: " "4371"
## [1] "s: " "4372"
## [1] "s: " "4373"
## [1] "s: " "4374"
## [1] "s: " "4375"
## [1] "s: " "4376"
## [1] "s: " "4377"
## [1] "s: " "4378"
## [1] "s: " "4379"
## [1] "s: " "4380"
## [1] "s: " "4381"
## [1] "s: " "4382"
## [1] "s: " "4383"
## [1] "s: " "4384"
## [1] "s: " "4385"
## [1] "s: " "4386"
## [1] "s: " "4387"
## [1] "s: " "4388"
## [1] "s: " "4389"
## [1] "s: " "4390"
## [1] "s: " "4391"
## [1] "s: " "4392"
## [1] "s: " "4393"
## [1] "s: " "4394"
## [1] "s: " "4395"
## [1] "s: " "4396"
## [1] "s: " "4397"
## [1] "s: " "4398"
## [1] "s: " "4399"
## [1] "s: " "4400"
## [1] "s: " "4401"
## [1] "s: " "4402"
## [1] "s: " "4403"
## [1] "s: " "4404"
## [1] "s: " "4405"
## [1] "s: " "4406"
## [1] "s: " "4407"
## [1] "s: " "4408"
## [1] "s: " "4409"
## [1] "s: " "4410"
## [1] "s: " "4411"
## [1] "s: " "4412"
## [1] "s: " "4413"
## [1] "s: " "4414"
## [1] "s: " "4415"
## [1] "s: " "4416"
## [1] "s: " "4417"
## [1] "s: " "4418"
## [1] "s: " "4419"
## [1] "s: " "4420"
## [1] "s: " "4421"
## [1] "s: " "4422"
## [1] "s: " "4423"
## [1] "s: " "4424"
## [1] "s: " "4425"
## [1] "s: " "4426"
## [1] "s: " "4427"
## [1] "s: " "4428"
## [1] "s: " "4429"
## [1] "s: " "4430"
## [1] "s: " "4431"
## [1] "s: " "4432"
## [1] "s: " "4433"
## [1] "s: " "4434"
## [1] "s: " "4435"
## [1] "s: " "4436"
## [1] "s: " "4437"
## [1] "s: " "4438"
## [1] "s: " "4439"
## [1] "s: " "4440"
## [1] "s: " "4441"
## [1] "s: " "4442"
## [1] "s: " "4443"
## [1] "s: " "4444"
## [1] "s: " "4445"
## [1] "s: " "4446"
## [1] "s: " "4447"
## [1] "s: " "4448"
## [1] "s: " "4449"
## [1] "s: " "4450"
## [1] "s: " "4451"
## [1] "s: " "4452"
## [1] "s: " "4453"
## [1] "s: " "4454"
## [1] "s: " "4455"
## [1] "s: " "4456"
## [1] "s: " "4457"
## [1] "s: " "4458"
## [1] "s: " "4459"
## [1] "s: " "4460"
## [1] "s: " "4461"
## [1] "s: " "4462"
## [1] "s: " "4463"
## [1] "s: " "4464"
## [1] "s: " "4465"
## [1] "s: " "4466"
## [1] "s: " "4467"
## [1] "s: " "4468"
## [1] "s: " "4469"
## [1] "s: " "4470"
## [1] "s: " "4471"
## [1] "s: " "4472"
## [1] "s: " "4473"
## [1] "s: " "4474"
## [1] "s: " "4475"
## [1] "s: " "4476"
## [1] "s: " "4477"
## [1] "s: " "4478"
## [1] "s: " "4479"
## [1] "s: " "4480"
## [1] "s: " "4481"
## [1] "s: " "4482"
## [1] "s: " "4483"
## [1] "s: " "4484"
## [1] "s: " "4485"
## [1] "s: " "4486"
## [1] "s: " "4487"
## [1] "s: " "4488"
## [1] "s: " "4489"
## [1] "s: " "4490"
## [1] "s: " "4491"
## [1] "s: " "4492"
## [1] "s: " "4493"
## [1] "s: " "4494"
## [1] "s: " "4495"
## [1] "s: " "4496"
## [1] "s: " "4497"
## [1] "s: " "4498"
## [1] "s: " "4499"
## [1] "s: " "4500"
## [1] "s: " "4501"
## [1] "s: " "4502"
## [1] "s: " "4503"
## [1] "s: " "4504"
## [1] "s: " "4505"
## [1] "s: " "4506"
## [1] "s: " "4507"
## [1] "s: " "4508"
## [1] "s: " "4509"
## [1] "s: " "4510"
## [1] "s: " "4511"
## [1] "s: " "4512"
## [1] "s: " "4513"
## [1] "s: " "4514"
## [1] "s: " "4515"
## [1] "s: " "4516"
## [1] "s: " "4517"
## [1] "s: " "4518"
## [1] "s: " "4519"
## [1] "s: " "4520"
## [1] "s: " "4521"
## [1] "s: " "4522"
## [1] "s: " "4523"
## [1] "s: " "4524"
## [1] "s: " "4525"
## [1] "s: " "4526"
## [1] "s: " "4527"
## [1] "s: " "4528"
## [1] "s: " "4529"
## [1] "s: " "4530"
## [1] "s: " "4531"
## [1] "s: " "4532"
## [1] "s: " "4533"
## [1] "s: " "4534"
## [1] "s: " "4535"
## [1] "s: " "4536"
## [1] "s: " "4537"
## [1] "s: " "4538"
## [1] "s: " "4539"
## [1] "s: " "4540"
## [1] "s: " "4541"
## [1] "s: " "4542"
## [1] "s: " "4543"
## [1] "s: " "4544"
## [1] "s: " "4545"
## [1] "s: " "4546"
## [1] "s: " "4547"
## [1] "s: " "4548"
## [1] "s: " "4549"
## [1] "s: " "4550"
## [1] "s: " "4551"
## [1] "s: " "4552"
## [1] "s: " "4553"
## [1] "s: " "4554"
## [1] "s: " "4555"
## [1] "s: " "4556"
## [1] "s: " "4557"
## [1] "s: " "4558"
## [1] "s: " "4559"
## [1] "s: " "4560"
## [1] "s: " "4561"
## [1] "s: " "4562"
## [1] "s: " "4563"
## [1] "s: " "4564"
## [1] "s: " "4565"
## [1] "s: " "4566"
## [1] "s: " "4567"
## [1] "s: " "4568"
## [1] "s: " "4569"
## [1] "s: " "4570"
## [1] "s: " "4571"
## [1] "s: " "4572"
## [1] "s: " "4573"
## [1] "s: " "4574"
## [1] "s: " "4575"
## [1] "s: " "4576"
## [1] "s: " "4577"
## [1] "s: " "4578"
## [1] "s: " "4579"
## [1] "s: " "4580"
## [1] "s: " "4581"
## [1] "s: " "4582"
## [1] "s: " "4583"
## [1] "s: " "4584"
## [1] "s: " "4585"
## [1] "s: " "4586"
## [1] "s: " "4587"
## [1] "s: " "4588"
## [1] "s: " "4589"
## [1] "s: " "4590"
## [1] "s: " "4591"
## [1] "s: " "4592"
## [1] "s: " "4593"
## [1] "s: " "4594"
## [1] "s: " "4595"
## [1] "s: " "4596"
## [1] "s: " "4597"
## [1] "s: " "4598"
## [1] "s: " "4599"
## [1] "s: " "4600"
## [1] "s: " "4601"
## [1] "s: " "4602"
## [1] "s: " "4603"
## [1] "s: " "4604"
## [1] "s: " "4605"
## [1] "s: " "4606"
## [1] "s: " "4607"
## [1] "s: " "4608"
## [1] "s: " "4609"
## [1] "s: " "4610"
## [1] "s: " "4611"
## [1] "s: " "4612"
## [1] "s: " "4613"
## [1] "s: " "4614"
## [1] "s: " "4615"
## [1] "s: " "4616"
## [1] "s: " "4617"
## [1] "s: " "4618"
## [1] "s: " "4619"
## [1] "s: " "4620"
## [1] "s: " "4621"
## [1] "s: " "4622"
## [1] "s: " "4623"
## [1] "s: " "4624"
## [1] "s: " "4625"
## [1] "s: " "4626"
## [1] "s: " "4627"
## [1] "s: " "4628"
## [1] "s: " "4629"
## [1] "s: " "4630"
## [1] "s: " "4631"
## [1] "s: " "4632"
## [1] "s: " "4633"
## [1] "s: " "4634"
## [1] "s: " "4635"
## [1] "s: " "4636"
## [1] "s: " "4637"
## [1] "s: " "4638"
## [1] "s: " "4639"
## [1] "s: " "4640"
## [1] "s: " "4641"
## [1] "s: " "4642"
## [1] "s: " "4643"
## [1] "s: " "4644"
## [1] "s: " "4645"
## [1] "s: " "4646"
## [1] "s: " "4647"
## [1] "s: " "4648"
## [1] "s: " "4649"
## [1] "s: " "4650"
## [1] "s: " "4651"
## [1] "s: " "4652"
## [1] "s: " "4653"
## [1] "s: " "4654"
## [1] "s: " "4655"
## [1] "s: " "4656"
## [1] "s: " "4657"
## [1] "s: " "4658"
## [1] "s: " "4659"
## [1] "s: " "4660"
## [1] "s: " "4661"
## [1] "s: " "4662"
## [1] "s: " "4663"
## [1] "s: " "4664"
## [1] "s: " "4665"
## [1] "s: " "4666"
## [1] "s: " "4667"
## [1] "s: " "4668"
## [1] "s: " "4669"
## [1] "s: " "4670"
## [1] "s: " "4671"
## [1] "s: " "4672"
## [1] "s: " "4673"
## [1] "s: " "4674"
## [1] "s: " "4675"
## [1] "s: " "4676"
## [1] "s: " "4677"
## [1] "s: " "4678"
## [1] "s: " "4679"
## [1] "s: " "4680"
## [1] "s: " "4681"
## [1] "s: " "4682"
## [1] "s: " "4683"
## [1] "s: " "4684"
## [1] "s: " "4685"
## [1] "s: " "4686"
## [1] "s: " "4687"
## [1] "s: " "4688"
## [1] "s: " "4689"
## [1] "s: " "4690"
## [1] "s: " "4691"
## [1] "s: " "4692"
## [1] "s: " "4693"
## [1] "s: " "4694"
## [1] "s: " "4695"
## [1] "s: " "4696"
## [1] "s: " "4697"
## [1] "s: " "4698"
## [1] "s: " "4699"
## [1] "s: " "4700"
## [1] "s: " "4701"
## [1] "s: " "4702"
## [1] "s: " "4703"
## [1] "s: " "4704"
## [1] "s: " "4705"
## [1] "s: " "4706"
## [1] "s: " "4707"
## [1] "s: " "4708"
## [1] "s: " "4709"
## [1] "s: " "4710"
## [1] "s: " "4711"
## [1] "s: " "4712"
## [1] "s: " "4713"
## [1] "s: " "4714"
## [1] "s: " "4715"
## [1] "s: " "4716"
## [1] "s: " "4717"
## [1] "s: " "4718"
## [1] "s: " "4719"
## [1] "s: " "4720"
## [1] "s: " "4721"
## [1] "s: " "4722"
## [1] "s: " "4723"
## [1] "s: " "4724"
## [1] "s: " "4725"
## [1] "s: " "4726"
## [1] "s: " "4727"
## [1] "s: " "4728"
## [1] "s: " "4729"
## [1] "s: " "4730"
## [1] "s: " "4731"
## [1] "s: " "4732"
## [1] "s: " "4733"
## [1] "s: " "4734"
## [1] "s: " "4735"
## [1] "s: " "4736"
## [1] "s: " "4737"
## [1] "s: " "4738"
## [1] "s: " "4739"
## [1] "s: " "4740"
## [1] "s: " "4741"
## [1] "s: " "4742"
## [1] "s: " "4743"
## [1] "s: " "4744"
## [1] "s: " "4745"
## [1] "s: " "4746"
## [1] "s: " "4747"
## [1] "s: " "4748"
## [1] "s: " "4749"
## [1] "s: " "4750"
## [1] "s: " "4751"
## [1] "s: " "4752"
## [1] "s: " "4753"
## [1] "s: " "4754"
## [1] "s: " "4755"
## [1] "s: " "4756"
## [1] "s: " "4757"
## [1] "s: " "4758"
## [1] "s: " "4759"
## [1] "s: " "4760"
## [1] "s: " "4761"
## [1] "s: " "4762"
## [1] "s: " "4763"
## [1] "s: " "4764"
## [1] "s: " "4765"
## [1] "s: " "4766"
## [1] "s: " "4767"
## [1] "s: " "4768"
## [1] "s: " "4769"
## [1] "s: " "4770"
## [1] "s: " "4771"
## [1] "s: " "4772"
## [1] "s: " "4773"
## [1] "s: " "4774"
## [1] "s: " "4775"
## [1] "s: " "4776"
## [1] "s: " "4777"
## [1] "s: " "4778"
## [1] "s: " "4779"
## [1] "s: " "4780"
## [1] "s: " "4781"
## [1] "s: " "4782"
## [1] "s: " "4783"
## [1] "s: " "4784"
## [1] "s: " "4785"
## [1] "s: " "4786"
## [1] "s: " "4787"
## [1] "s: " "4788"
## [1] "s: " "4789"
## [1] "s: " "4790"
## [1] "s: " "4791"
## [1] "s: " "4792"
## [1] "s: " "4793"
## [1] "s: " "4794"
## [1] "s: " "4795"
## [1] "s: " "4796"
## [1] "s: " "4797"
## [1] "s: " "4798"
## [1] "s: " "4799"
## [1] "s: " "4800"
## [1] "s: " "4801"
## [1] "s: " "4802"
## [1] "s: " "4803"
## [1] "s: " "4804"
## [1] "s: " "4805"
## [1] "s: " "4806"
## [1] "s: " "4807"
## [1] "s: " "4808"
## [1] "s: " "4809"
## [1] "s: " "4810"
## [1] "s: " "4811"
## [1] "s: " "4812"
## [1] "s: " "4813"
## [1] "s: " "4814"
## [1] "s: " "4815"
## [1] "s: " "4816"
## [1] "s: " "4817"
## [1] "s: " "4818"
## [1] "s: " "4819"
## [1] "s: " "4820"
## [1] "s: " "4821"
## [1] "s: " "4822"
## [1] "s: " "4823"
## [1] "s: " "4824"
## [1] "s: " "4825"
## [1] "s: " "4826"
## [1] "s: " "4827"
## [1] "s: " "4828"
## [1] "s: " "4829"
## [1] "s: " "4830"
## [1] "s: " "4831"
## [1] "s: " "4832"
## [1] "s: " "4833"
## [1] "s: " "4834"
## [1] "s: " "4835"
## [1] "s: " "4836"
## [1] "s: " "4837"
## [1] "s: " "4838"
## [1] "s: " "4839"
## [1] "s: " "4840"
## [1] "s: " "4841"
## [1] "s: " "4842"
## [1] "s: " "4843"
## [1] "s: " "4844"
## [1] "s: " "4845"
## [1] "s: " "4846"
## [1] "s: " "4847"
## [1] "s: " "4848"
## [1] "s: " "4849"
## [1] "s: " "4850"
## [1] "s: " "4851"
## [1] "s: " "4852"
## [1] "s: " "4853"
## [1] "s: " "4854"
## [1] "s: " "4855"
## [1] "s: " "4856"
## [1] "s: " "4857"
## [1] "s: " "4858"
## [1] "s: " "4859"
## [1] "s: " "4860"
## [1] "s: " "4861"
## [1] "s: " "4862"
## [1] "s: " "4863"
## [1] "s: " "4864"
## [1] "s: " "4865"
## [1] "s: " "4866"
## [1] "s: " "4867"
## [1] "s: " "4868"
## [1] "s: " "4869"
## [1] "s: " "4870"
## [1] "s: " "4871"
## [1] "s: " "4872"
## [1] "s: " "4873"
## [1] "s: " "4874"
## [1] "s: " "4875"
## [1] "s: " "4876"
## [1] "s: " "4877"
## [1] "s: " "4878"
## [1] "s: " "4879"
## [1] "s: " "4880"
## [1] "s: " "4881"
## [1] "s: " "4882"
## [1] "s: " "4883"
## [1] "s: " "4884"
## [1] "s: " "4885"
## [1] "s: " "4886"
## [1] "s: " "4887"
## [1] "s: " "4888"
## [1] "s: " "4889"
## [1] "s: " "4890"
## [1] "s: " "4891"
## [1] "s: " "4892"
## [1] "s: " "4893"
## [1] "s: " "4894"
## [1] "s: " "4895"
## [1] "s: " "4896"
## [1] "s: " "4897"
## [1] "s: " "4898"
## [1] "s: " "4899"
## [1] "s: " "4900"
## [1] "s: " "4901"
## [1] "s: " "4902"
## [1] "s: " "4903"
## [1] "s: " "4904"
## [1] "s: " "4905"
## [1] "s: " "4906"
## [1] "s: " "4907"
## [1] "s: " "4908"
## [1] "s: " "4909"
## [1] "s: " "4910"
## [1] "s: " "4911"
## [1] "s: " "4912"
## [1] "s: " "4913"
## [1] "s: " "4914"
## [1] "s: " "4915"
## [1] "s: " "4916"
## [1] "s: " "4917"
## [1] "s: " "4918"
## [1] "s: " "4919"
## [1] "s: " "4920"
## [1] "s: " "4921"
## [1] "s: " "4922"
## [1] "s: " "4923"
## [1] "s: " "4924"
## [1] "s: " "4925"
## [1] "s: " "4926"
## [1] "s: " "4927"
## [1] "s: " "4928"
## [1] "s: " "4929"
## [1] "s: " "4930"
## [1] "s: " "4931"
## [1] "s: " "4932"
## [1] "s: " "4933"
## [1] "s: " "4934"
## [1] "s: " "4935"
## [1] "s: " "4936"
## [1] "s: " "4937"
## [1] "s: " "4938"
## [1] "s: " "4939"
## [1] "s: " "4940"
## [1] "s: " "4941"
## [1] "s: " "4942"
## [1] "s: " "4943"
## [1] "s: " "4944"
## [1] "s: " "4945"
## [1] "s: " "4946"
## [1] "s: " "4947"
## [1] "s: " "4948"
## [1] "s: " "4949"
## [1] "s: " "4950"
## [1] "s: " "4951"
## [1] "s: " "4952"
## [1] "s: " "4953"
## [1] "s: " "4954"
## [1] "s: " "4955"
## [1] "s: " "4956"
## [1] "s: " "4957"
## [1] "s: " "4958"
## [1] "s: " "4959"
## [1] "s: " "4960"
## [1] "s: " "4961"
## [1] "s: " "4962"
## [1] "s: " "4963"
## [1] "s: " "4964"
## [1] "s: " "4965"
## [1] "s: " "4966"
## [1] "s: " "4967"
## [1] "s: " "4968"
## [1] "s: " "4969"
## [1] "s: " "4970"
## [1] "s: " "4971"
## [1] "s: " "4972"
## [1] "s: " "4973"
## [1] "s: " "4974"
## [1] "s: " "4975"
## [1] "s: " "4976"
## [1] "s: " "4977"
## [1] "s: " "4978"
## [1] "s: " "4979"
## [1] "s: " "4980"
## [1] "s: " "4981"
## [1] "s: " "4982"
## [1] "s: " "4983"
## [1] "s: " "4984"
## [1] "s: " "4985"
## [1] "s: " "4986"
## [1] "s: " "4987"
## [1] "s: " "4988"
## [1] "s: " "4989"
## [1] "s: " "4990"
## [1] "s: " "4991"
## [1] "s: " "4992"
## [1] "s: " "4993"
## [1] "s: " "4994"
## [1] "s: " "4995"
## [1] "s: " "4996"
## [1] "s: " "4997"
## [1] "s: " "4998"
## [1] "s: " "4999"
## [1] "s: " "5000"
pred_mean_healthy <- unlist(map(MEAN, mean))
pred_df_healthy = tibble(
time = t,
pred_rate = pred_mean_healthy,
emp_rate = y
)
pred_df_healthy %>%
ggplot() +
geom_line(aes(x = time, y = emp_rate), col = "blue") +
geom_line(aes(x = time, y = pred_rate), col = "orange")
par(mfrow = c(3, 2))
for(j in sample(12, 3)){
plot(ALPHA[[j]], type = 'l', ylab = paste("Alpha", j), xlab = "Iter", col = "purple",
main = paste("Acceptance ratio =", round(length(unique(ALPHA[[j]])) / S, 3)))
acf(ALPHA[[j]], ylab = paste("Alpha", j), main = paste("Series Alpha", j))
}
par(mfrow = c(3, 2))
for(j in sample(17, 3)){
plot(DELTA[[j]], type = 'l', ylab = paste("Delta", j), xlab = "Iter", col = "purple",
main = paste("Acceptance ratio =", round(length(unique(DELTA[[j]])) / S, 3)))
acf(DELTA[[j]], ylab = paste("Delta", j), main = paste("Series Delta", j))
}
par(mfrow = c(1,2))
plot(DC, type = 'l', ylab = "DC", xlab = "Iter", col = "purple",
main = paste("Acceptance ratio =", round(length(unique(DC)) / S, 3)))
acf(DC, ylab = "DC", main = "Series DC")
par(mfrow = c(1,2))
plot(TAO, type = 'l', ylab = "Tao", xlab = "Iter", col = "purple",
main = paste("Acceptance ratio =", round(length(unique(TAO)) / S, 3)))
acf(TAO, ylab = "Tao", main = "Series Tao")
alpha_names <- map_chr(1:12, function(x) str_c("alpha_", x, sep = ""))
alpha_df <- map_dfc(ALPHA, function(x) x)
colnames(alpha_df) <- alpha_names
alpha_df <- alpha_df %>%
pivot_longer(1:12, names_to = "index", values_to = "alpha_1") %>%
mutate(index = factor(index, levels = alpha_names))
colnames(alpha_df)[2] <- "value"
alpha_df %>%
filter(index %in% alpha_names[1:6]) %>%
ggplot() +
geom_density(aes(x = value, fill = index), alpha = 0.7)
alpha_df %>%
filter(index %in% alpha_names[7:12]) %>%
ggplot() +
geom_density(aes(x = value, fill = index), alpha = 0.7)
delta_names <- map_chr(1:17, function(x) str_c("delta_", x, sep = ""))
delta_df <- map_dfc(DELTA, function(x) x)
colnames(delta_df) <- delta_names
delta_df <- delta_df %>%
pivot_longer(1:17, names_to = "index", values_to = "delta_1") %>%
mutate(index = factor(index, levels = delta_names))
colnames(delta_df)[2] <- "value"
delta_df %>%
filter(index %in% delta_names[1:6]) %>%
ggplot() +
geom_density(aes(x = value, fill = index), alpha = 0.7)
delta_df %>%
filter(index %in% delta_names[7:12]) %>%
ggplot() +
geom_density(aes(x = value, fill = index), alpha = 0.7)
delta_df %>%
filter(index %in% delta_names[13:17]) %>%
ggplot() +
geom_density(aes(x = value, fill = index), alpha = 0.7)
for(j in sample(length(MEAN), 6)){
plot_2d <- ggplot() +
geom_point(aes(x = MEAN[[j]], y = TAO)) +
geom_density_2d(aes(x = MEAN[[j]], y = TAO))
print(plot_2d)
}
samples_arr <- read_csv("samples.csv", skip = 1)
## Parsed with column specification:
## cols(
## `'seconds'` = col_double(),
## `'mV'` = col_double()
## )
colnames(samples_arr) <- c("time", "signal")
y <- samples_arr$signal[76:662] #signal 0.208 max(samples$signal[0:250]) which(samples$signal[0:250] == 0.58)
t <- samples_arr$time[76:662] #signal 0.705 max(samples$signal[0:700]) which(samples$signal[0:700] == 0.705)
S <- 5000
burn <- 5000
beta_curr <- 1
alpha_curr <- rep(1, 12)
delta_curr<- seq(0, max(t), length.out = 17)
end_time <- max(t)
tao_curr <- 1
mean_curr <- rep(mean(y), length(y)) #mean_curr is a vector cause mean_star of y is a vector so I think r only lets me run it if it is a vector
dc_curr <- 0.5
miu_curr <- mean_curr - dc_curr
#storage
DC <- NULL
ALPHA <- rep(list(NULL), 12)
BETA <- NULL
DELTA <- rep(list(NULL), 17)
TAO <- NULL
MEAN <- rep(list(NULL), length(y))
for(s in 1:S) {
#print(c("s: ", s))
#1 draw of Beta
beta_lst <- beta_one_samp(beta_curr, BETA, s)
beta_curr <- beta_lst[[1]]
BETA <- beta_lst[[2]]
#print(c("beta: ", beta_curr))
#1 draw of Alphas (1-12)
for (k in 1:12) {
alpha_lst <- alpha_one_samp(k, alpha_curr[k], ALPHA[[k]], s)
alpha_curr[k] <- alpha_lst[[1]]
ALPHA[[k]] <- alpha_lst[[2]]
}
#print(c("ALPHA: ", alpha_curr))
#1 draw of Deltas (1-17)
delta_lst <- delta_one_samp(delta_curr[1], 0, delta_curr[2], DELTA[[1]], s)
delta_curr[1] <- delta_lst[[1]]
DELTA[[1]] <- delta_lst[[2]]
for (k in 2:16) {
delta_lst <- delta_one_samp(delta_curr[k], delta_curr[k-1], delta_curr[k+1], DELTA[[k]], s)
delta_curr[k] <- delta_lst[[1]]
DELTA[[k]] <- delta_lst[[2]]
}
delta_lst <- delta_one_samp(delta_curr[17], delta_curr[[16]], end_time, DELTA[[17]], s)
delta_curr[17] <- delta_lst[[1]]
DELTA[[17]] <- delta_lst[[2]]
#print(c("DELTA: ", delta_curr))
#1 draw of Tao
tao_lst <- tao_one_samp_cond(tao_curr, TAO, s)
tao_curr <- tao_lst[[1]]
TAO <- tao_lst[[2]]
#print(c("TAO: ", tao_curr))
#1 draw of DC
dc_lst <- dc_one_samp(dc_curr, DC, s)
dc_curr <- dc_lst[[1]]
DC <- dc_lst[[2]]
#print(c("DC: ", dc_curr))
mean_curr <- dc_curr + gen_miu_star(t, alpha_curr = alpha_curr, beta_curr = beta_curr, delta_curr = delta_curr)
for (k in 1:length(y)) {
MEAN[[k]] <- c(MEAN[[k]], mean_curr[[k]])
}
#print(c("MEAN CURR: ", mean_curr))
}
pred_mean_arr <- unlist(map(MEAN, mean))
pred_df_arr = tibble(
time = t,
pred_rate = pred_mean_arr,
emp_rate = y
)
pred_df_arr %>%
ggplot() +
geom_line(aes(x = time, y = emp_rate), col = "blue") +
geom_line(aes(x = time, y = pred_rate), col = "orange")
par(mfrow = c(3, 2))
for(j in sample(12, 3)){
plot(ALPHA[[j]], type = 'l', ylab = paste("Alpha", j), xlab = "Iter", col = "purple",
main = paste("Acceptance ratio =", round(length(unique(ALPHA[[j]])) / S, 3)))
acf(ALPHA[[j]], ylab = paste("Alpha", j), main = paste("Series Alpha", j))
}
par(mfrow = c(3, 2))
for(j in sample(17, 3)){
plot(DELTA[[j]], type = 'l', ylab = paste("Delta", j), xlab = "Iter", col = "purple",
main = paste("Acceptance ratio =", round(length(unique(DELTA[[j]])) / S, 3)))
acf(DELTA[[j]], ylab = paste("Delta", j), main = paste("Series Delta", j))
}
par(mfrow = c(1,2))
plot(DC, type = 'l', ylab = "DC", xlab = "Iter", col = "purple",
main = paste("Acceptance ratio =", round(length(unique(DC)) / S, 3)))
acf(DC, ylab = "DC", main = "Series DC")
par(mfrow = c(1,2))
plot(TAO, type = 'l', ylab = "Tao", xlab = "Iter", col = "purple",
main = paste("Acceptance ratio =", round(length(unique(TAO)) / S, 3)))
acf(TAO, ylab = "Tao", main = "Series Tao")
alpha_names <- map_chr(1:12, function(x) str_c("alpha_", x, sep = ""))
alpha_df <- map_dfc(ALPHA, function(x) x)
colnames(alpha_df) <- alpha_names
alpha_df <- alpha_df %>%
pivot_longer(1:12, names_to = "index", values_to = "alpha_1") %>%
mutate(index = factor(index, levels = alpha_names))
colnames(alpha_df)[2] <- "value"
alpha_df %>%
filter(index %in% alpha_names[1:6]) %>%
ggplot() +
geom_density(aes(x = value, fill = index), alpha = 0.7)
alpha_df %>%
filter(index %in% alpha_names[7:12]) %>%
ggplot() +
geom_density(aes(x = value, fill = index), alpha = 0.7)
delta_names <- map_chr(1:17, function(x) str_c("delta_", x, sep = ""))
delta_df <- map_dfc(DELTA, function(x) x)
colnames(delta_df) <- delta_names
delta_df <- delta_df %>%
pivot_longer(1:17, names_to = "index", values_to = "delta_1") %>%
mutate(index = factor(index, levels = delta_names))
colnames(delta_df)[2] <- "value"
delta_df %>%
filter(index %in% delta_names[1:6]) %>%
ggplot() +
geom_density(aes(x = value, fill = index), alpha = 0.7)
delta_df %>%
filter(index %in% delta_names[7:12]) %>%
ggplot() +
geom_density(aes(x = value, fill = index), alpha = 0.7)
delta_df %>%
filter(index %in% delta_names[13:17]) %>%
ggplot() +
geom_density(aes(x = value, fill = index), alpha = 0.7)
for(j in sample(length(MEAN), 6)){
plot_2d <- ggplot() +
geom_point(aes(x = MEAN[[j]], y = TAO)) +
geom_density_2d(aes(x = MEAN[[j]], y = TAO))
print(plot_2d)
}
